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ABSTRACT 

We investigate the use of a wide variety of spectroscopic measurements to determine distances to low-redshift Type la supernovae 
(SN la) in the Hubble flow observed through the CfA Supernova Program. We consider linear models for predicting distances to SN 
la using light-curve width and color parameters (determined using the SALT2 light-curve fitter) and a spectroscopic indicator, and 
evaluate the resulting Hubble diagram scatter using a cross-validation procedure. We confirm the ability of spectral flux ratios alone 
at maximum light to reduce the scatter of Hubble residuals by ~ 10% [weighted rms, or WRMS = 0.189 ± 0.026 mag for the flux 
ratio '/?(6630/4400)] with respect to the standard combination of light-curve width and color, for which WRMS = 0.204 ± 0.029 mag. 
When used in combination with the SALT2 color parameter, the color-corrected flux ratio ^?'^(6420/5290) at maximum light leads 
to an even lower scatter (WRMS = 0.175 ± 0.025 mag), although the improvement has low statistical significance (< 2cr) given the 
size of our sample (26 SN la). We highlight the importance of an accurate relative flux calibration and the failure of this method for 
highly-reddened objects. Comparison with synthetic spectra from 2D delayed-detonation explosion models shows that the correlation 
of 7^(6630/4400) with SN la absolute magnitudes can be largely attributed to intrinsic color variations and not to reddening by 
dust in the host galaxy. We consider flux ratios at other ages, as well as the use of pairs of flux ratios, revealing the presence of 
small-scale intrinsic spectroscopic variations in the iron-group-dominated absorption features around ~ 4300 A and ~ 4800 A. The 
best flux ratio overall is the color-corrected 7?'' (46 10/4260) at f = -2.5 d from maximum light, which leads to ~ 30% lower scatter 
(WRMS = 0.143 ± 0.020 mag) with respect to the standard combination of light-curve width and color, at ~ 2cr significance. We 
examine other spectroscopic indicators related to line-profile morphology (absorption velocity, pseudo-equivalent width etc.), but 
none appear to lead to a significant improvement over the standard light-curve width and color parameters. We discuss the use of 
spectra in measuring more precise distances to SN la and the implications for future surveys which seek to determine the properties 
of dark energy. 

Key words, supernovae: general — cosmology: observations 



1. Introduction mological constant (for which w = -1) within about 10% (Astier 
\ I et al. 2006; Wood-Vasey et al. 2007). Constraints on the varia- 
J> Precise distances to Type la supernovae (SN la) formed the cor- tion of w with redshift come from from high-redshift observa- 
nerstone of the discovery of cosmic acceleration (Riess et al. tions with the Hubble Space Telescope (Riess et al. 2004, 2007). 
1998; Perlmutteret al. 1999). These measurements use the shape Present-day limits are weak, but future work with large, care- 
of supernova light curves and their colors to tell which super- fully calibrated samples from the ground (Pan-STARRS, Dark 
novae are bright and which are intrinsically dim (Phillips 1993; Energy Survey, LSST) and from space (Euclid, WFIRST) will 
Riess et al. 1996; Prieto et al. 2006; Jha et al. 2007; Guy et al. contribute to distinguishing the nature of dark energy (Albrecht 
2007; Conley et al. 2008; Mandel et al. 2009). In this paper we et al. 2009). In designing the follow-up observations for these en- 
explore the suggestion of Bailey et al. (2009) that spectra can terprises, it is worth knowing whether spectra will be useful only 
contribute to improved distance measurements. We apply statis- for classification and precise redshifts, or whether the spectra of 
tical tests to a subset of the ~ 250 SN la for which we have the supernovae themselves can be used to improve the precision 
good light curves and spectra based on the ongoing program of of the distances. The way we explore this is to analyze the CfA 
supernova observations at the Harvard-Smithsonian Center for sample, using the difference between the distance derived from 
Astrophysics (CfA; Matheson et al. 2008; Hicken et al. 2009a). Hubble expansion with the distance predicted from our various 
It is important to construct the best possible distance indi- "^o^^^ls- difference is the Hubble residual, which we use as 
cators to extract the maximum cosmological information from ^ ^^^^^^^ of the power of a particular model to predict the su- 
supernova surveys. The present state-of-the-art gives distances P^™"^^ '^'^t^'^'^e- described below, we explore models that 
to well-observed individual objects with uncertainties of order '^""^^ine quantitative information from the spectrum with mfor- 
10%, so that samples of nearby (Hicken et al. 2009b) and distant "^^^1°" °" l^S^t curve shape and color. 

SN la (ESSENCE, Miknaitis et al. 2007; SNLS, Astier et al. Spectroscopic information is fundamental to the success of 

2006) can be combined to constrain the equation-of-state for employing SN la as distance indicators in large surveys. Cleanly 

dark energy, noted w. The first results show that for a flat uni- separating Type la supernovae from core-collapse events like 

verse with constant w, the dark energy is compatible with a cos- SN lb and SN Ic improves the purity of the sample. More di- 
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rectly, Nugent et al. (1995) showed that some easily-measured 
line ratios in SN la spectra are correlated with the luminos- 
ity. Measurements of line velocities (and gradients thereof), 
strengths, and widths and their relation to supernova luminos- 
ity have been explored recently by several authors (Benetti et al. 
2005; Blondin et al. 2006; Bongard et al. 2006; Hachinger et al. 
2006; Bronder et al. 2008). Likewise, Matheson et al. (2008) re- 
vealed spectroscopic variability amongst SN la of similar lumi- 
nosity. But the first application of spectroscopic clues to improve 
distance estimates has come from Bailey et al. (2009). Using 
spectra of 58 SN la from the Nearby Supernova Factory, they 
showed that the ratio of fluxes in selected wavelength bins (flux 
ratios) could reduce the scatter of Hubble residuals by ~ 20% 
compared to the usual combination of light-curve width and 
color parameters (cr = 0.128±0.012mag cf. 0.161 ±0.015mag). 
By using a flux ratio measured on a de-reddened spectrum in 
combination with a color parameter they found a further ~ 5% 
improvement (cr = 0.119 + 0.01 1 mag). We have sought first to 
see if we can reproduce their results using the CfA data set, and 
then to test additional ideas about ways to use spectra to improve 
the estimates of supernova distances. 

In practice, the standardization of SN la magnitudes involves 
a term related to the width of the light curve and a correction 
due to color. While some methods attempt to separate intrin- 
sic color variations from reddening by dust in the host galaxy 
(e.g., MLCS2k2; Jha et al. 2007) others use a single parameter 
for both eff'ects (e.g. SALT2; Guy et al. 2007), exploiting the 
degeneracy between the two: underluminous SN la are also in- 
trinsically redder than overluminous SN la (e.g. Tripp 1998). We 
adopt the latter approach in this paper, to match the method used 
by Bailey et al. (2009). An active area of research involves the 
use of SN la spectra to provide independent or complementary 
information on SN la luminosities that would help improve their 
use as distance indicators. 

We consider models for predicting distances to SN la of the 
form: 



We present the CfA data set in § 3. In § 4 we study the flux 
ratios of Bailey et al. (2009), while in § 5 we consider other 
spectroscopic indicators. We discuss the use of SN la spectra for 
distance measurements in § 6 and conclude in § 7. 



2. Methodology 

2. 1 . Light-curve fitting 

We use the S ALT2 hght-curve fitter of Guy et al. (2007) to deter- 
mine the width and color parameters for each SN la in our sam- 
ple. A model relating distance, apparent magnitudes, and hnear 
dependencies of the absolute magnitude is: 



fi — Mb — M + axi — 13c + yS, 



(2) 



where (xi , c) are the S ALT2 light curve width and color parame- 
ters, and S is some spectroscopic indicator. The rest-frame peak 
apparent fi-band magnitude is m^, also obtained from the SALT2 
fit to a supernova's Ught curve. The distance modulus predicted 
from the light curve and spectral indicators is yU, and the constant 
M is a reference absolute magnitude. The distance modulus es- 
timated from the redshift is yu(z) = 25 + 5 logio[DL(z)Mpc~'] 
under a fixed cosmology, where is the luminosity distance. 

We use the exact same SALT2 options as Guy et al. (2007) 
to fit the SN la light curves in our sample, and only trust the re- 
sult when the following conditions are met: reduced xl ^ 2; at 
least one B-band point before +5 d from B-band maximum, and 
one after h-10 d; at least 5 B- and V-band points in the age range 
-15 <t< -1-60 d; finally, we impose a cut on the SALT2 xi pa- 
rameter, namely -3 < xi < 2. This last condition is equivalent 
to considering SN la in the range 0.8 < Ami5(B) < 1.7 (i.e. sub- 
luminous 1991bg-like SN la are excluded). We examined aU the 
Ught-curve fits by eye to ensure they were satisfactory given this 
set of conditions. Approximately 170 of the ~ 250 SN la with 
Ught curves from the CfA SN program pass these requirements. 



H - Mb - M + (a X width) - (fix color) + (y x spec), (1) 2.2. Training 



where is the apparent rest- frame B-band magnitude at peak, 
M is a reference absolute magnitude, "width" and "color" are the 
usual light-curve parameters, and "spec" is some spectroscopic 
indicator; {a,/3,j) are fitting constants. We study the following 
five models: 

1. only a spectroscopic indicator is used [i.e. ia,P) = (0, 0)], 

2. both a spectroscopic indicator and a light-curve width pa- 
rameter are used, but no color parameter (i.e. /3 - 0), 

3. both a spectroscopic indicator and a color parameter are 
used, but no light-curve width parameter (i.e. a = 0), 

4. a spectroscopic indicator is used in addition to the light- 
curve width and color parameters. 

5. both light-curve width and color parameters are used, but no 
spectroscopic indicator (i.e. y = 0). We refer to this as the 
"standard" model. 

We refer to the set of light-curve parameters and spectroscopic 
indicators in a given model as the "predictors" for that model, 
as is common practice in the field of statistics. We can evaluate 
the use of including a spectroscopic indicator (models 1-4) by 
comparing the resulting scatter of Hubble diagram residuals with 
that from the standard model (No. 5). 

The paper is organized as follows: in § 2 we present our light- 
curve fitting and training method, as well as a cross-validation 
procedure to evaluate the impact of each spectroscopic indicator. 



For estimating the coefficients of the model (training), we 
use a custom version of the luminosity distance fitter 
simple.cosf itter' (A. Conley 2009, private communication) 
based on the Minuit function minimization package (James & 
Roos 1975). This code minimizes the foUowing expression with 
respect to the parameters (a,jS, y, M): 



2 _ V ^'^B,s - >npred,siZs;xl,Cs,Ss;a,/3,y,M)f 

S= 1 

_ Y iMimg, x{,c„ Sj, a,fi, y, M) - n{zs)f 

.v=l * 



(3) 



where mB,s is the rest- frame peak apparent 5-band magnitude of 
the 5* SN la, and mpred.s is the predicted peak apparent B-band 
magnitude, given by: 

OTpred,* = 5 logio Dl{Zs) - CUXl^s + PCs - ySs + M, (4) 

where Dl is the luminosity distance at redshift z, for a 
given cosmological model described by the standard parameters 
(w, Qm, Qa)- Since our analysis only includes objects at low red- 
shifts (z < 0.06), we do not solve for these parameters and sim- 
ply assume a flat, cosmological constant-dominated model with 



' http://qold.astro.utoronto.ca/conley/simple_cosfitter 
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(w,Q„,nA) = (-1,0.27,0.73). The M term is a collection of 
constants including the reference M. 

The variance cr^ that appears in the denominator of Eq. 3 in- 
cludes an error on the corrected magnitude (Eq. 2), using the 
estimation error covariance of the light-curve parameters and 
spectroscopic indicators, a variance due to peculiar velocities 
[Cpec,.! - (vpec/cz.s)(5/ In 10), whcrc we take the rms peculiar 
velocity Vpec - 300 km s '], and an intrinsic dispersion of SN la 
magnitudes: 

(5) 



2 2 2 



where CTmt is adjusted iteratively until xl ~ ^ [typically (Tint ^ 
0.2 mag for the standard (xi , c) model] . For any particular model, 
the intrinsic variance cr^ , accounts for deviations in magnitude 

int 

in the Hubble diagram beyond that explained by measurement 
error or random peculiar velocities, and hence represents a floor 
to how accurately the model can predict distances. To limit the 
impact of the peculiar velocity error we restrict our analysis to 
SN la at redshifts z > 0.015 (o-p^c < 0.15 mag). Of the 170 SN la 
with satisfactory SALT2 fits, 114 are at redshifts greater than 
0.015. 



2.3. Cross-validation 

We consider several models described by Eq. 2 that use differ- 
ent subsets of the predictors ix\,c,S). If we train a model on 
the data of all the SN in the sample to estimate the coefficients 
ia,0,f,M), we can evaluate the fit of the model by comput- 
ing the training error, e.g. the mean squared distance modulus 
residual, yu(OTg, x^, Cs, Ss', a,P, J, M) - fiiZs), over all SN 5 in the 
training set. 

For finite samples, the average Hubble diagram residual of 
the training set SN is an optimistic estimate of the abiUty of 
the statistical model, Eq. 2, to make accurate predictions given 
the supernova observables. This is because it uses the supernova 
data twice: first for estimating the model parameters (training), 
and second in evaluating the residual error. Hence, the training 
set residuals underestimate the prediction error, which is the ex- 
pected error in estimating the distance of a SN that was not orig- 
inally in the finite training set. We refer to these data as "out- 
of-sample". Furthermore, with a fixed, finite, and noisy training 
data set, it is always possible to reduce the residual, or train- 
ing, error of the fit by introducing more predictors to the model. 
However, this may lead to over-fitting, in which apparently sig- 
nificant predictors are found in noisy data, even though in reality 
there was no trend. These relationships are sensitive to the finite 
training set and would not generalize to out-of-sample cases. To 
evaluate predictive performance and guard against over-fitting 
with a statistical model based on finite data, we should estimate 
the prediction error for out-of-sample cases. To do so, we use 
a cross-validation (CV) procedure to evaluate the impact of us- 
ing a spectroscopic indicator S, alone and in conjunction with 
standard light curve parameters, on the accuracy of distance pre- 
dictions in the Hubble diagram. 

Cross-validation seeks to estimate prediction error and to test 
the sensitivity of the trained statistical model to the data set by 
partitioning the full data set into smaller subsets. One subset is 
held out for testing predictions of the model, while its comple- 
ment is used to train the model. This process is repeated over 
partitions of the full data set. This method avoids using the same 
data simultaneously for training the model and for estimating 
its prediction error. Cross-validation was used before for statisti- 
cal modeling of SN la by Mandel et al. (2009), who applied the 
.632 bootstrap method to evaluate distance predictions for SN la 



using near infrared light curves. A careful implementation of a 
cross-validation method is particularly important for small sam- 
ples, as is the case in this paper (e.g. 26 SN la at maximum light; 
see § 4.3). 

In this paper, the cross-validation method we use is known 
as ^T-fold CV. The idea is to divide our SN la sample into K 
subsets, train a given model on K - I subsets, and validate it 
on the remaining subset. This procedure is repeated K times, at 
which point all SN la have been part of a validation set once. 
Typical choices of are 5 or 10 (e.g., Hastie et al. 2009). The 
case K = N, where A'^ is the number of SN la in our sample, is 
known as "leave-one-out" CV. In this case, each SN la in turn is 
used as a validation set, and the training is repeated A'^ times on 
N-ISN la. 

In practice, we run /iT-fold CV as follows: 

1. the sample of SN la is randomly divided into K subsets 
of equal size (when A'^ is not a multiple of K, the number of 
SN la between any two subsets differs by at most one). 

2. Looping over each K fold: 

2a. all the SN la in the A"* subset are removed from the sam- 
ple: they form the validation set. The remaining SN la 
define the training set. 

2b. the objects in the training set are then used to determine 
the best-fit values for the parameters (ax, 0k, Jk, Mk) in 
Eq. 3, as well as the intrinsic dispersion (Tint in Eq. 5. 

2c. using this set of parameters we predict the magnitudes of 
the SN la in the validation set (indexed j): 



5\og^(^DL(zj)-aKXij+l3KCj-yKSj+MK- (6) 



The Hubble residual, or error, of the predicted distance 
modulus is then 

/S.Hj = mBj-m^rcdj = K^B, x{,Cj, Sf, a,p, y, M) -fiizj). 

(7) 

3. When the magnitude or distance of each SN la has been pre- 
dicted once using the above scheme, we analyze the predic- 
tion errors (§2.4). When doing so, we check that the set of 
best-fit {uk^Pk, 7k, Mk) are consistent amongst all training 
sets. 

For all the spectroscopic indicators we consider in this paper, 
we run /T-fold CV with K = 2,5,10, and A^ to make sure our re- 
sults are not sensitive to the exact choice of K (the impact on the 
weighted rms of prediction Hubble residuals is < 0.002 mag). 
Moreover, we run each K-fo\d CV 10 times to check the out- 
come is insensitive to how the starting SN la sample is divided 
into K subsets (the impact on the weighted rms of prediction 
Hubble residuals is < 0.003 mag). In what follows we report our 
results based onK = 10. 



2.4. Comparing modei predictions 

For each model, which we label by its predictors, e.g. (xi , c, S), 
cross-validation gives us a set of prediction errors {A//, j for each 
SN s. To summarize the total dispersion of predictions, we com- 
puted the weighted mean squared error, 



WRMS^ 



s=l 



(8) 



s=l 



the square root of which is the weighted rms. We weight the con- 
tribution from each SN by the inverse of its expected total vari- 
ance (the precision) Ws = cr^^. We prefer to use the rms of the 
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prediction residuals rather than the sample standard deviation, 
since the former measures the average squared deviation of the 
distance prediction from the Hubble distance ju(z), whereas the 
latter measures the average squared deviation of prediction er- 
rors from the mean prediction error. Note that the mean squared 
error is equal to the sample variance plus the square of the mean 
error. Thus, the mean squared error will be larger than the sam- 
ple variance if the mean error, or bias, is significant, but the two 
statistics will be the same if it is not. Since the mean prediction 
error is not guaranteed to be zero, we use the WRMS statistic to 
assess the total dispersion of distance prediction errors. We also 
estimate the sampling variance of this statistic (see Appendix A). 

The WRMS measures the total dispersion in the Hubble 
diagram. However, we expect that some of that scatter is due 
to random peculiar velocities [influencing fi(z) with variance 
CTpgj,], and some due to measurement error (cr^,). Using the cross- 
vahdated distance errors, we also estimate how precisely we can 
expect a particular model to predict the distance to a SN la when 
these other sources of error are negligible. We call this vari- 
ance estimate the rms intrinsic prediction error, a property of 
the model itself, and label it o-^^^^- Intuitively, this is the result of 
subtracting from the total dispersion the expected contributions 
of peculiar velocities and measurement uncertainties. It is simi- 
lar to the intrinsic variance criai discussed in § 2.2, in that it rep- 
resents a floor to how accurately the model can predict distances. 
It is not strictly equivalent, however, since CTmt is adjusted dur- 
ing the training process so that;(fj « 1, while cr^tA is estimated 
using the cross-validated distance modulus prediction errors. In 
Appendix B, we describe a maximum likelihood estimate for 
cTpred and its standard error from the set of distance predictions. 

We are also interested in the intrinsic covciricince of the 
distance prediction errors generated by two different models. 
Imagine that peculiar velocities and measurement error were 
negligible, and model P and model Q predict distances to the 
same set of SN la. We calculate the prediction errors, {Ajt/^, AyUp) 
from each model. There is a positive intrinsic covariance if A//f 
tends to be positive when A//^ is positive, and a negative intrinsic 
covariance if they tend to make errors in opposite directions. The 
intrinsic correlation is important because it suggests how useful 
it would be to combine the distance predictions of two models. If 
two models tend to make prediction errors in the same direction 
(positive correlation), then the combined model is not likely to 
do much better than the most accurate of the two original mod- 
els. However, if two models tend to make prediction errors that 
are wrong in different ways (zero or negative correlation), then 
we expect to see a gain from averaging the two models. 

Even if two models make prediction errors that are intrin- 
sically uncorrelated, random pecuUar velocities will tend to in- 
duce a positive correlation in the realized errors {A/if,AyUp) if 
the methods are used on the same set of SN. This is because the 
unknown pecuhar velocity for a given SN is the same regardless 
of the model we use to generate its distance prediction. Hence, 
the expected contribution of random peculiar velocities to the 
sample covariance of predictions must be removed to estimate 
the intrinsic covariance between two models. In Appendix B, 
we describe a maximum hkehhood estimator for the intrinsic 
covariance and its standard error using the set of distance pre- 
dictions. 

We use the maximum likelihood estimation method to es- 
timate the intrinsic prediction error and intrinsic covariance of 
each model compared to the reference model {x\,c) that uses 
only light curve information. 



3. Spectroscopic data 

We have used a large spectroscopic data set obtained through the 
CfA Supernova Program. Since 1994, we have obtained ~ 2400 
optical spectra of ~ 450 low-redshift (z < 0.05) SN la with the 
1.5 m Tillinghast telescope at FLWO using the FAST spectro- 
graph (Fabricant et al. 1998). Several spectra were published in 
studies of specific supernovae (e.g., SN 1998bu; Jha et al. 1999), 
while 432 spectra of 32 SN la have recently been published by 
Matheson et al. (2008). We also have complementary multi-band 
optical photometry for a subset of ~ 250 SN la (Riess et al. 
1999; Jha et al. 2006; Hicken et al. 2009a), as well as NIR JHKs 
photometry for the brighter ones (Wood-Vasey et al. 2008). All 
pubhshed data are available via the CfA Supernova Archive^. 

All the spectra were obtained with the same telescope and 
instrument, and reduced in a consistent manner (see Matheson 
et al. 2008 for details). The uniformity of this data set is unique 
and enables an accurate estimate of our measurement errors. 



4. Spectral flux ratios 

4.1. Measurements 

Bailey et al. (2009) introduced a new spectroscopic indicator, 
calculated as the ratio of fluxes in two wavelength regions of a 
SN la spectrum binned on a logarithmic wavelength scale. This 
ratio, noted 'R(Ax/Ay) = F(Ax)/F(Ay) [Ax and Ay being the rest- 
frame wavelength coordinates in A of a given bin center], is mea- 
sured on a de-redshifted spectrum corrected for Galactic redden- 
ing using the CardelU et al. (1989) extinction law with Ry = 3.1 
in combination with the dust maps of Schlegel et al. (1998). A 
color-corrected version of this flux ratio, noted '??''(/(x//If), is 
measured on a spectrum additionally corrected for the SALT2 
color parameter using the color law of Guy et al. (2007). Figure 1 
illustrates both measurements. 

We use the same binning as Bailey et al. (2009), namely 134 
bins equally spaced in In A between 3500 A and 8500 A (rest 
frame), although most of the CfA spectra used here do not ex- 
tend beyond ~ 7100 A (see § 4.3). The resulting ~ 2000 km s'^ 
bin size is significantly less than the typical width of a SN la 
feature (~ 10000 km s"'). The error on includes a flux error 
(from the corresponding variance spectrum), an error due to the 
relative flux calibration accuracy (see § 4.2), and an error due 
to the SALT2 color precision. When there are several spectra of 
a given SN la within Af = 2.5 d of the age we consider (see 
§ 4.3 for spectra at maximum light; § 4.4 for spectra at other 
ages), we use the error-weighted mean and standard deviation of 
all flux ratios as our measurement and error, respectively. Bailey 
et al. (2009) also chose At - 2.5 d in their analysis, and we find 
that increasing Af worsens the results while decreasing it leads 
to too small a sample. 

Bailey et al. (2009) cross-checked the results for their 
best single flux ratio '??(6420/4430) using the the sample of 
SN la spectra published by Matheson et al. (2008) [and avail- 
able through the CfA SN Archive]. We checked the validity 
of our flux ratio measurements by comparing the values of 
'??(6420/4430) in the Matheson et al. (2008) sample wifli fliose 
reported in Table 2 of Bailey et al. (2009). In all cases, our 
measurements agree well within the Icr errors. This also holds 
for SN 1998bu, accidentally removed from the Matheson et al. 
(2008) sample by Bailey et al. (2009) [H. Fakhouri 2010, private 
communication]. We note that we were unable to cross-check 

^ http://www.cfa.harvard.edu/supernova/SNarchive.html 
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2.5 - 



I ' ■ ■ I I ' ' ' I I ' ' ■ I ' ' I 
Ordinal spectrum: R 
(conccted for MW extinction) 

tag 




j I I I I I I I I ' i M I r I I ' i 1 ' i I I I I I I I I I I I I I I r 

De-reddeaed spectrmn: R° 
(corrected for SALT2 color) 




|A(B-V)| [nmg] 



4000 4500 5000 5500 6000 6500 
Log-binned Rest Wavelength [A] 

Fig. 1. Illustration of the flux ratio measurement. The upper panel 
shows the input spectrum (de-redshifted and corrected for Galactic red- 
dening; here SN 1999gd around maximum light), binned on a log- 
wavelength scale. The gray vertical lines represent the fluxes in charac- 
teristic wavelength bins mentioned throughout the text. The lower panel 
shows the same spectrum corrected for SALT2 color, which is used to 
measure the color-corrected flux ratios K^. 



the flux ratio measurements of Bailey et al. (2009) in a similar 
fashion, since none of their 58 SN la spectra are publicly avail- 
able. 



4.2. Impact of relative flux calibration and SALT2 color 

When the Ax and Ay wavelength bins have a large separation 
(> lOOOA), 'R{Ax/Ay) is essentially a color measurement. We 
therefore expect flux ratios to be sensitive to the relative flux 
calibration accuracy of the spectra. Fig. 2 shows the relation be- 
tween uncorrected Hubble residuals [i.e. - M - fj.(z)] and our 
most highly-ranked flux ratio ??(6630/4400) at maximum light 
(see § 4.3). There is one data point per SN la, color-coded ac- 
cording to the absolute diff'erence in Z? - V color at maximum 
light derived from the spectrum and that derived from the pho- 
tometry, noted \A(B - V)\, which we use as a proxy for rela- 
tive flux calibration accuracy. The bulk of the sample defines a 
highly correlated relation (dashed line), with several outliers all 
having \A{B - V)\ > 0.1 mag. We therefore restrict our analysis 
to SN la with spectra that have a relative flux calibration better 
than 0.1 mag. 

Bailey et al. (2009) noted that the highly-reddened 
SN 1999cl was a large outlier in their analysis, and attributed 
this to the non-standard nature of the extinction towards this SN 
(Rv ~ 1.5; Krisciunas et al. 2006). To explore the effects of red- 
dening, in Fig. 3 (left), we show the relation between uncorrected 
Hubble residual and '/?(6630/4400) at maximum light, for SN la 
at redshifts z > 0.005 that satisfy our requirement on the relative 
flux calibration accuracy. Using this lower redshift bound has 
the effect of including several highly-reddened SN la (including 
SN 1999cl; see Fig. 4), which are otherwise excluded based on 
the redshift cut we use elsewhere this paper (z > 0.015). For 
SN la with c < 0.5, ^^(6630/4400) is highly correlated with un- 



0.00 



0.07 0.13 020 0.26 0.33 0.39 




0.4 0.6 0.8 

Hux ratio R(6630j'44O0) 

Fig. 2. Uncorrected Hubble residual vs. flux ratio 7?(6630/4400) at 
maximum light, color-coded according to the absolute difference in 
B-V color derived from the spectrum and that derived from the photom- 
etry, noted |A(B - V)\. The dashed line is a linear fit to the SN la with 
|A(B - V)\ < 0.1 mag. The highly-reddened SN 2006br is not shown 
here. 
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Fig. 4. Histogram of the S ALT2 color parameter (c) for SN la at z > 
0.015 (open) and 0.005 <z< 0.015 {hatched). Bins that include SN la 
with c > 0.5 are labeled. 



corrected Hubble residuals, but those with red colors (c > 0.5) 
tend to deviate significantly from this relation (dashed line; this 
is not the case for SN 1995E, for which c ^ 0.9), the two 
largest outliers corresponding to the reddest SN la (SN 1999cl 
and SN 2006X). Both are subject to high extinction by non- 
standard dust in their respective host galaxies (Ay ~ 2 mag for 
Rv ~ 1.5; Krisciunas et al. 2006; Wang et al. 2008) and dis- 
play time-variable Na I D absorption, whose circumstellar or in- 
terstellar origin is still debated (Patat et al. 2007; Blondin et al. 
2009). The reddening curves in Fig. 3 (dotted lines) seem to cor- 
roborate the fact that the nonlinear increase of flux ratios at high 
values of the S ALT2 color parameter is mainly due to reddening 
by dust with low Ry. Nonetheless, SN 1999cl still stands out in 
this respect as it would require a value of Ry < 0.5 inconsistent 
with that found by Krisciunas et al. (2006). Moreover, while we 
obtain consistent Ry estimates for SN 2006X using other flux 
ratios, this is not the case for SN 2006br, for which some flux 
ratios are consistent with Rv - 3.1. 
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Fig. 3. Left: Uncorrected Hubble residual vs. flux ratio ';?(6630/4400) at maximum light for SN la at z > 0.005 with |A(S - V)| < 0.1 mag, 
color-coded according to the SALT2 color parameter, c. Points corresponding to SN la with c > 0.5 are labeled. The dashed line is a linear fit to 
the SN la with c < 0.5. The dotted lines are reddening curves for different values of Ry, normalized to the smallest '/?(6630/4400) value. Right: 
Color-corrected Hubble residual vs. color-corrected flux ratio '??''(6420/5290) at maximum light. 



The right panel of Fig. 3 shows the relation between color- 
corrected Hubble residual [i.e. ms - M -/3c-fi{z)] and our most 
highly-ranked color-corrected flux ratio ^'(6420/5290) at max- 
imum light (see § 4.3) for the same sample. SN la with a SALT2 
color c > 0.5 are again outliers. As noted by Bailey et al. 2009, 
this shows that a single color parameter cannot encompass the 
variety of SN la intrinsic colors and extinction by non-standard 
dust. We therefore impose a cut on SALT2 color in our analy- 
sis, only considering SN la with c < 0.5. Four of the five SNe 
with c > 0.5 in Fig. 3 are rejected anyway based on our redshift 
cut. The remaining one, SN 2006br, is then rejected based on our 
color cut. 



4.3. Results on Maximum-light Spectra 
4.3.1. Selecting the best flux ratios 

After selecting SN la that satisfy both requirements on relative 
flux calibration accuracy and SALT2 color parameter, we are left 
with 26 SN la at z > 0.015 with spectra within At - 2.5 d from 
maximum light (see Table 1, where we also present selected 
flux ratio measurements). The spectra show no sign of signifi- 
cant contamination by host-galaxy light, which can also bias the 
flux ratio measurements. We make no cut based on the signal- 
to-noise ratio (S/N) of our spectra, as they are generally well in 
excess of 100 per log- wavelength bin. We only consider flux ra- 
tios for wavelength bins represented in all the spectra. This leads 
to 98 bins between ~ 3690 A and ~ 7060 A, i.e. 9506 indepen- 
dent flux ratios. 

We run the /T-fold cross-validation procedure outlined in 
§ 2.3, and consider the five models for estimating distances to 
SN la described in § 1: 





- niB - M + yR 


(9) 


fl 


- niB - M + ax] + yR 


(10) 




= mB-M-pc + yl^ 


(11) 




- niB - M + axi - /3c + yj^ 


(12) 




— Mb — M + ax\ — /3c. 


(13) 



When no color correction is involved (Eqs. 9-10), we use the 
uncorrected flux ratio K. When a color correction is involved 
(the -ySc term in Eqs. 11-12), we use the color-corrected version 
of the flux ratio 'R'^. Using K in combination with color, or 'R'^ 
alone or in combination with xi , severely degrades the predictive 
power of the model, so we do not report results using (c, ■??); 'R'^ 
alone; or {xu'R'^). 

We rank the flux ratios in each case based on the intrinsic 
prediction error (o-pi-ed; see § 2.4), but note that ranking based 
on the weighted rms of prediction Hubble residuals makes al- 
most no difference. The results for the top five flux ratios are 
displayed in Table 2. We also report the best-fit y, the weighted 
rms of prediction Hubble residuals (WRMS), the intrinsic cor- 
relation of residuals with those found using the standard {xi , c) 
predictors (noted Pa,,c; see § 2.4), and the difference in intrinsic 
prediction error with respect to the standard (xi,c) model, noted 
A^i c. Since we compute the error on A^, ^ (see Appendix B), we 
also report the significance of this difference with respect to the 
standard (xi, c) predictors. This is a direct measure of whether a 
particular model predicts more accurate distances to SN la when 
compared to the standard approach, and if so how significant 
is the improvement. Fig. 5 shows the resulting Hubble diagram 
residuals vs. redshift for the best flux ratio in each of the four 
models given by Eqs. 9-12, and using the standard (jci,c) pre- 
dictors. 

All the flux ratios listed in Table 2 lead to an improvement 
over the standard (xi,c) correction (i.e. Av,,c < 0), as found by 
Bailey et al. (2009), but the significance is low: < Icr for R only; 
<K lo- for (xi,R.y, ~ 1.5cr for (cH'') and (xi,c,R''). This is in 
part due to the small number of SN la in our sample. Note that 
Pai,c > 0.5 in all cases, i.e. the models that include a flux ratio 
tend to make prediction errors in the same direction as (^i , c), 
and we do not expect to gain much by combining these models. 

Using best single flux ratio '^(6630/4400) by itself reduces 
the weighted rms of prediction residuals (as well as the intrinsic 
prediction error, (Xpied) by < 10% when compared with (xi,c) 
[WRMS = 0.189 + 0.026 mag cf. 0.204 + 0.029 mag], although 
as noted above the significance of the difference in intrinsic pre- 
diction error is negligible (Av|,e = -0.018 +0.025 mag, or 0.7cr). 
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Table 1. SN la sample for flux ratio measurements at maximum light 



SN zcMB ms 5 c ?;(6630/4400) ?y(6420/5290) ?y(5690/5360) ?y(5160/5290) ^^'CSegO/SSSO) 



1998V 


0.0170 


15.085 (0.020) 


-0.435 (0.161) 


0.031 (0.015) 


0.330 (0.005) 


0.744 (0.006) 


1.045 (0.004) 


1.485 (0.006) 


0.933 (0.004) 


1998dx 


0.0539 


17.536 (0.037) 


-1.890(0.457) 


-0.051 (0.027) 


0.365 (0.018) 


0.949 (0.027) 


1.086 (0.024) 


1.445 (0.033) 


0.946 (0.018) 


1998eg 


0.0237 


16.096 (0.016) 


-0.956 (0.366) 


0.048 (0.019) 


0.378 (0.023) 


0.920 (0.022) 


1.070(0.011) 


1.478 (0.013) 


0.995 (0.007) 


1999aa 


0.0152 


14.698 (0.009) 


0.896 (0.073) 


-0.019 (0.009) 


0.286 (0.027) 


0.704 (0.025) 


1.054 (0.012) 


1.448 (0.011) 


0.994 (0.006) 


1999CC 


0.0316 


16.760 (0.010) 


-1.891 (0.175) 


0.057 (0.012) 


0.406 (0.013) 


0.994 (0.017) 


1.001 (0.015) 


1.513(0.019) 


0.850 (0.011) 


1999ek 


0.0176 


15.587 (0.009) 


-1.075 (0.127) 


0.163 (0.010) 


0.478 (0.039) 


0.936 (0.036) 


1.091 (0.016) 


1.461 (0.014) 


0.952 (0.007) 


1999gd 


0.0191 


16.940 (0.022) 


-1.210(0.193) 


0.455 (0.022) 


0.714 (0.013) 


0.929 (0.016) 


1.081 (0.024) 


1.455 (0.020) 


1.028 (0.015) 


2000dk 


0.0165 


15.347 (0.021) 


-2.658 (0.301) 


0.055 (0.022) 


0.423 (0.017) 


0.933 (0.017) 


1.038 (0.009) 


1.425 (0.010) 


1.006 (0.006) 


2000fa 


0.0218 


15.883 (0.023) 


0.311 (0.127) 


0.100(0.018) 


0.410 (0.046) 


0.829 (0.042) 


1.082(0.020) 


1.459 (0.017) 


0.902 (0.009) 


2001 eh 


0.0363 


16.575 (0.018) 


1.457 (0.222) 


0.020 (0.017) 


0.312(0.047) 


0.799 (0.043) 


1.040(0.020) 


1.451 (0.019) 


0.931 (0.010) 


2002ck 


0.0302 


16.303 (0.048) 


-0.183 (0.147) 


-0.017 (0.023) 


0.348 (0.046) 


0.752 (0.042) 


1.030(0.020) 


1.463 (0.017) 


0.945 (0.010) 


2002hd 


0.0360 


16.738 (0.038) 


-0.748 (0.456) 


0.100 (0.022) 


0.403 (0.027) 


0.748 (0.027) 


0.955 (0.018) 


1.378 (0.020) 


0.899 (0.013) 


2002hu 


0.0359 


16.587 (0.012) 


0.052 (0.143) 


-0.052 (0.012) 


0.293 (0.027) 


0.768 (0.029) 


1.048 (0.021) 


1.509(0.029) 


0.984 (0.017) 


2002jy 


0.0187 


15.702 (0.019) 


0.660 (0.212) 


0.013 (0.015) 


0.305 (0.015) 


0.795 (0.015) 


1.073 (0.010) 


1.566(0.012) 


0.963 (0.007) 


2002kf 


0.0195 


15.654 (0.033) 


-1.493 (0.189) 


0.009 (0.023) 


0.361 (0.035) 


0.953 (0.035) 


1.132(0.024) 


1.550 (0.024) 


0.969 (0.014) 


2003U 


0.0279 


16.471 (0.046) 


-2.536 (0.558) 


0.033 (0.035) 


0.373 (0.009) 


0.891 (0.016) 


1.048(0.016) 


1.445 (0.025) 


0.965 (0.014) 


2003ch 


0.0256 


16.659 (0.022) 


-1.655 (0.297) 


0.012 (0.019) 


0.402 (0.025) 


1.024(0.024) 


1.218(0.010) 


1.570(0.009) 


1.043 (0.005) 


2003it 


0.0240 


16.342 (0.028) 


-1.815 (0.359) 


0.084 (0.029) 


0.432(0.011) 


0.918(0.011) 


1.087(0.007) 


1.456(0.009) 


0.950 (0.005) 


2003iv 


0.0335 


16.961 (0.026) 


-2.473 (0.486) 


-0.031 (0.028) 


0.420 (0.005) 


1.092(0.011) 


1.081 (0.011) 


1.467 (0.015) 


0.940 (0.009) 


2004as 


0.0321 


16.956 (0.018) 


-0.017 (0.206) 


0.128 (0.016) 


0.415 (0.006) 


0.838 (0.008) 


1.065(0.009) 


1.573 (0.011) 


0.975 (0.007) 


2005ki 


0.0208 


15.551 (0.029) 


-2.123 (0.153) 


-0.059 (0.026) 


0.398 (0.032) 


0.975 (0.030) 


1.091 (0.015) 


1.441 (0.015) 


0.957 (0.009) 


2006ax 


0.0180 


15.010 (0.010) 


-0.062 (0.062) 


-0.049 (0.009) 


0.304 (0.019) 


0.851 (0.018) 


1.082(0.008) 


1.558 (0.009) 


1.009(0.005) 


2006gj 


0.0277 


17.668 (0.033) 


-2.073 (0.280) 


0.409 (0.023) 


0.675 (0.008) 


0.988 (0.012) 


1.011 (0.011) 


1.377 (0.014) 


0.930 (0.009) 


2006sr 


0.0232 


16.126(0.017) 


-1.754 (0.220) 


0.060 (0.015) 


0.422 (0.019) 


0.939 (0.018) 


1.031 (0.010) 


1.474 (0.010) 


0.913 (0.006) 


2007ca 


0.0151 


15.933 (0.013) 


0.289 (0.122) 


0.305 (0.012) 


0.547 (0.028) 


0.895 (0.026) 


1.097(0.012) 


1.569 (0.012) 


0.986 (0.006) 


2008bf 


0.0257 


15.703 (0.010) 


0.097 (0.095) 


0.031 (0.010) 


0.315 (0.038) 


0.806 (0.035) 


1.047(0.016) 


1.526 (0.013) 


0.924 (0.007) 



Notes. Spectra within 2.5 d from maximum light for these SN la are available via the Cf A Supernova Archive. 



Table 2. Top 5 flux ratios at maximum Ught from 10-fold CV on 26 SN la 



Rank 


Ax 


Ay 


r 


WRMS 




Pxi,c 






1 


6630 


4400 


-4.37 ± 0.09 


0.189 ±0.026 


0.163 ±0.030 


0.80 ± 0.09 


-0.018 ±0.025 (O.70-) 


2 


6630 


4430 


-4.44 + 0.12 


0.191 ±0.027 


0.166 ±0.030 


0.60 


±0.15 


-0.015 ±0.033 (0.5O-) 


3 


6630 


4670 


-5.16 + 0.11 


0.197 ±0.027 


0.171 ±0.031 


0.57 


±0.17 


-0.007 ± 0.032 (0.2O-) 


4 


6900 


4460 


-5.71 ±0.12 


0.196 ±0.027 


0.171 ±0.031 


0.54 


±0.17 


-0.009 ± 0.036 (0.2cr) 


5 


6420 


4430 


-3.40 ±0.10 


0.196 ±0.028 


0.173 ±0.031 


0.56 


±0.16 


-0.010 ± 0.033 (0.3O-) 




1 


6630 


4400 


-4.51+0.15 


0.201 ± 0.028 


0.176 + 0.032 


0.74 + 0.10 


-0.006 + 0.028 (O.20-) 


2 


6630 


4430 


-4.57 ± 0.20 


0.204 ± 0.029 


0.180 ±0.032 


0.55 ±0.17 


-0.001 ± 0.033 (O.Ocr) 


3 


6900 


4460 


-5.81 ±0.16 


0.206 ± 0.029 


0.182 ±0.032 


0.52 


±0.17 


0.002 ±0.032 (O.lcr) 


4 


6900 


4370 


-4.46 ± 0.04 


0.206 ± 0.029 


0.182 ±0.032 


0.69 


±0.13 


0.001 ± 0.030 (0.0a-) 


5 


6990 


4370 


-4.95 + 0.05 


0.207 + 0.029 


0.183 + 0.032 


0.61+0.15 


0.003 + 0.035(0.10-) 




1 


6420 


5290 


-1.75 ±0.10 


0.175 ±0.025 


0.148 ±0.029 


0.80 


±0.09 


-0.032 ±0.023 (1.4a-) 


2 


6630 


4890 


-2.19 ±0.14 


0.181 ±0.024 


0.148 ±0.031 


0.81 


±0.08 


-0.023 ± 0.022 (1.0a-) 


3 


4890 


6630 


0.73 + 0.06 


0.182 ±0.024 


0.148 + 0.031 


0.83 + 0.08 


-0.023 + 0.023 (1.0a-) 


4 


4890 


6810 


0.60 + 0.05 


0.180 + 0.024 


0.148 + 0.030 


0.73 + 0.11 


-0.026 ±0.026 (1.0a-) 


5 


6540 


4890 


-1.84 ±0.12 


0.179 ±0.024 


0.149 ±0.030 


0.88 


±0.05 


-0.026 ±0.020 (1.3a-) 




1 


5690 


5360 


-2.78 ± 0.20 


0.164 ±0.023 


0.134 ±0.028 


0.69 


±0.13 


-0.044 ±0.028 (1.6cr) 


2 


5360 


5690 


3.23 ± 0.25 


0.167 ±0.023 


0.137 ±0.028 


0.69 


±0.13 


-0.042 ±0.028 (1.5cr) 


3 


5660 


5290 


-1.95 ±0.24 


0.173 ±0.024 


0.142 ±0.029 


0.72 


±0.11 


-0.038 ±0.029 (1.30-) 


4 


5690 


5290 


-1.67 ±0.21 


0.171 ±0.024 


0.142 ±0.028 


0.75 


±0.11 


-0.037 ±0.022 (1.70-) 


5 


5290 


5660 


2.75 ± 0.33 


0.174 ±0.024 


0.144 ±0.029 


0.73 


±0.11 


-0.036 ± 0.028 (1.30-) 



ixi,c) 



0.204 + 0.029 0.181 ±0.032 



Notes. WRMS is the error-weighted rms of prediction Hubble residuals (in magnitudes; see Appendix A); o^prej is the intrinsic prediction error (in 
magnitudes); Pxi,c is the intrinsic correlation in prediction error with the standard (xi , c) predictors; last, A^^^c is the difference in intrinsic prediction 
error with respect to {xi,c) [see Appendix B]. 



Using "R in combination with xi leads to no improvement 
over using alone (although this is not reported by Bailey et al. 
2009, it is consistent with their findings; S. Bailey 2009, pri- 
vate communication), and even leads to systematically worse 
results. Our best single flux ratio ??(6630/4400) yields a dif- 
ference in intrinsic prediction error with respect to (xi,c) of 
^xi,c = -0.018 ± 0.025 mag, when used on its own, while it 



yields A^^^^ = -0.006 ± 0.028 mag when combined with xi. 
These differences are statistically indistinguishable from one an- 
other given the size of the error on A i, _c, but they are systematic 
regardless of the flux ratio we consider. 

This seems counter-intuitive, as one might expect that in- 
cluding an additional predictor would result in more accurate 
distance predictions. However, this is not necessarily the case 
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Fig. 5. Hubble diagram residuals for the highest-ranked flux ratios at 
maximum light. From top to bottom: prediction residuals using "R only; 
CR, xi); ('??'', c); ('R'^, xj, c); and using the standard SALT2 fit parameters 
(xi,c). In each case we indicate the weighted rms of prediction Hubble 
residuals {gray highlighted region). 



under cross-validation. The reason is that xi by itself is a poor 
predictor of Hubble residuals, and one does not gain anything 
by combining it with '??(6630/4400). This is not surprising, as 
the relation between light-curve width and luminosity is only 
valid if the SN la are corrected for color or extinction by dust 
beforehand. In fact, !??(663 0/4400) by itself accounts for most of 
the variation in Hubble residuals. When we cross-validate, the 
extra coefficient a will tend to fit some noise in a given train- 
ing set, and this relation will not generahze to the validation set. 
This results in an increase in prediction error because the added 
information is not useful. We see from Table 2 that adding xi 
affects the best-fit value for y[y = -4.5 1 ± 0.15 cf. -4.37 ± 0.09 
for ??(6630/4400) only]; moreover, we obtain a < when us- 
ing [xi,??(6630/4400)] where a » 0.15 when using (xi,c), 
which again shows that a is fitting noise when ^(6630/4400) is 
combined with xi . This illustrates the advantage of using cross- 
validation in guarding against over-fitting noise as more param- 
eters and potential predictors are added. 

Figure 6 {upper panel) shows why ??(6630/4400) alone is a 
good predictor of Hubble residuals. Its strong correlation with 
SALT2 color (Pearson correlation coefficient r - 0.92) shows 
that this ratio is essentially a color measurement. The correlation 
with xi is less pronounced (r = —0.38), but this is largely due to 
a small number of outliers: removing the three largest outliers 
results in a Pearson correlation coefficient r = -0.65. The flux 
ratio '/?(6630/4400) by itself is thus as useful a predictor as xi 
and c combined. 
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Fig. 6. Correlation between the highest-ranked (??, Rf^) at maximum 
light and the SALT2 fit parameters (xi, c). 



The relation between ??(663 0/4400) and xi is not linear, 
but it is certainly true that SN la with higher xi (i.e. broader 
light curves) tend to have lower '??(6630/4400) [the same is true 
for !/?(6420/4430), the highest-ranked flux ratio by Bailey et al. 
2009]. Since the width of the Ughtcurve is a parameter intrinsic 
to each SN la (although its measurement can be subtly affected 
by host-galaxy reddening; see Phillips et al. 1999), the correla- 
tion between xi and '??(6630/4400) shows that the color variation 
measured by ??(6630/4400) is intrinsic in part. This is consis- 
tent with the so-called "brighter-bluer" relation of Tripp (1998): 
overluminous SN la are intrinsically bluer than underluminous 
SN la (see also Riess et al. 1996). 

Using a color-corrected flux ratio 'R.'^ in combination with 
color results in even lower Hubble residual scatter when com- 
pared with the single flux ratio case. Our best flux ratio in 
this case, '^''(6420/5290), reduces the weighted rms of predic- 
tion residuals by ~ 15% with respect to (xi,c) [WRMS = 
0.175 ± 0.025 mag cf. 0.204 ± 0.029mag], and the intrinsic 
prediction error by ~ 20% [cTpred = 0.148 ± 0.029 mag cf. 
0.181 ± 0.032 mag]. Again, the significance of this difference is 
only 1.4cr (A^,,,. = -0.032 + 0.023 mag). We see from Fig. 6 
{middle panel) that (6420/5290) is strongly anti-correlated 
with xi {r = -0.78), and that dereddening the spectra using 
the SALT2 color law is effective in removing any dependence 
of ??"(6420/5290) on color, as expected. 

One would naively think that combining our best color- 
corrected ratio 9?*^ (6420/5290) with (xi,c) would lead to an 
even further improvement, but this is not the case. In fact, 
!R^(6420/5290) ranks 298* when we consider the set of pre- 
dictors (xi,c,!R^). This is due to the strong anti-correlation of 
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'R'^ (6420/5290) with xi. Adding xi as an extra predictor when 
'^'■(6420/5290) akeady includes this information means a will 
tend to fit noise in a given training set, as was the case for the 
set of (xi , K) predictors when compared with !R-only. Indeed, the 
best-fit value for a for [xi, c, ^^"^ (6420/5290)] is again consistent 
with 0. 

Nonetheless, several color-corrected flux ratios do result in a 
further reduced scatter when combined with (xi , c), although the 
wavelength baseline for these ratios is much smaller (< 400 A) 
and the wavelength bins forming the ratios are all concentrated in 
the region of the 811/1/15454,5640 doublet. Our highest-ranked 
flux ratio in this case, ??'^(5690/5360), reduces the weighted 
rms of prediction residuals by ~ 20% with respect to (xi,c) 
[WRMS = 0.164 ± 0.023 mag cf. 0.204 + 0.029 mag], and the 
intrinsic prediction error by ~ 25% [cTpred = 0.134 ± 0.028 cf. 
0.181 + 0.032mag]. Again, the significance of this difference is 
only 1.60- (A.t,,c = -0.044 + 0.028 mag). We see from Fig. 6 
(lower panel) that this ratio is not correlated with xi (r = -0.08) 
or c (r = 0.11), and thus constitutes a useful additional predictor 
of distances to SN la. 
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Fig. 8. Weighted rms of prediction Hubble residuals v.s-. absolute 
Pearson cross-correlation coefficient for all flux ratios at maximum 
light, excluding (top) and including (bottom) the highly-reddened 
SN2006br. 



4.3.2. Two-dimensional maps of all flux ratios 

The results for all 9506 flux ratios are displayed in Fig. 7. The 
four rows correspond to the four models for estimating SN la 
distances that include a flux ratio (Eqs. 9-12). The left column is 
color-coded according to the weighted rms of prediction Hubble 
residuals (flux ratios that result in WRMS > 0.324 mag are given 
the color corresponding to WRMS = 0.324 mag), while the right 
colunm is color-coded according to the absolute Pearson correla- 
tion coeflicient of the correction terms with uncorrected Hubble 
residuals [e.g. for the set of predictors (c, 'R^), the correlation of 
(-j8c -I- YRf^) with uncorrected residuals]. 

Only a very restricted number of wavelength bins lead to 
a low WRMS of prediction Hubble residuals when a flux ratio 
"R. is used by itself (Fig. 7; upper left), namely Ax <: 6300 A 
and Ay * 4400 A (4 of flie 5 best flux ratios in Table 2 for K- 
only have Ay « 4400 A). This is in stark contrast with the large 
number of flux ratios with absolute Pearson correlation coeffi- 
cients \r\ > 0.8 (Fig. 7; upper right). In general, a flux ratio with 
a higher correlation coefficient will result in a Hubble diagram 
with less scatter, but this is not systematically the case, and the 
relation between the two is certainly not linear. For Pearson cor- 
relation coefficients \r\ > 0.8, the standard deviation of Hubble 
residuals can vary by up to 0. 1 mag at any given |r| (Fig. 8, top 
panel). This is because the cross-correlation coefficient does not 
take into account errors on K or on the Hubble residual, and is 
biased by outliers and reddened SN la. The lower panel of Fig. 8 
shows the impact of including the highly-reddened SN 2006br: 
at any given \r\, the resulting weighted rms of prediction Hubble 
residuals is 30-60% higher. Moreover, many flux ratios with high 
correlation coefficients (\r\ > 0.8) result in Hubble diagrams with 
excessively large scatter (WRMS > 1 mag). This is counter- 
intuitive, since the resulting scatter in these cases appears to be 
larger than when no predictors at all are used to determine dis- 
tances to SN la (in which case WRMS « 0.5 mag). The rea- 
son is that we consider the scatter under cross-validation, as op- 
posed to fitting all the SN la at the same time. In these aberrant 
cases, the trained model is sensitive to the inclusion or exclu- 
sion of some outUer in the training set, and this leads to large 
errors when the outlier is in the validation set. Last, including 
this SN leads to correlations with \r\ > 0.95, where there are 
none otherwise. Fig. 8 thus justifies our excluding SN 2006br 



from the sample (already excluded based on our cut on SALT2 
color; see § 4.2), and illustrates the advantage of selecting flux 
ratios based directly on the weighted rms of prediction Hubble 
diagram residuals, rather than on cross-correlation coefficients. 
As already mentioned in § 4.3.1, using cross-validated predic- 
tion errors to select the best flux ratios guards us against overfit- 
ting a small sample: in the naive approach that consists in fitting 
the entire SN la sample at once, adding more predictors always 
leads to a lower scatter in Hubble residuals (this is known as 
"resubstitution"; see, e.g., Mandel et al. 2009). 

When the SALT2 color parameter is used in combination 
with a color-corrected flux ratio 'R'^, there are again restricted 
wavelength regions that lead to a low weighted rms of predic- 
tion Hubble residuals (Fig. 7; third row left) [4 of the 5 best 
flux ratios in Table 2 for (c, 'R'^) predictors involve wavelength 
bins at ~ 4900 A and ~ 6500 A]. The SALT2 color parameter c 
does not attempt to distinguish between reddening by dust and 
intrinsic color variations. Dereddening the spectra using this pa- 
rameter corrects for both effects regardless of their relative im- 
portance. However, since the SALT2 color law is very similar 
to file CardeUi et al. (1989) extinction law with Ry = 3.1 and 
E(B - y) = 0.1 mag (Guy et al. 2007, their Fig. 3), one generally 
assumes that the color correction removes the bulk of reddening 
by dust, and the remaining variations in the SED are primarily 
intrinsic to the supernova. If this is so, it is intriguing that the 
best flux ratios for the !??-only and (c, 'R'^) models share similar 
wavelength bins. The recent survey of 2D SN la models from 
Kasen et al. (2009) suggests that a significant part of the color 
variation measured by the ^(6630/4400) is indeed intrinsic (see 
§ 4.3.4). 

The second row of Fig. 7 confirms that using the xi parame- 
ter in combination with a flux ratio results in a slight degradation 
in the weighted rms of prediction residuals, while the correla- 
tions with uncorrected Hubble residuals are degraded with re- 
spect to cases where R is used by itself. Last, the bottom row of 
Fig. 7 is a visual demonstration that (xi,c,'R^) fares better than 
(c, 'R'^) overall, although the best color-corrected flux ratios do 
not perform significantly better. We see from the right panel that 
the correlations of (axi - /3c + jR!^) with uncorrected residuals 
all have absolute Pearson correlation coefficients \r\ < 0.5. The 
two regions at Axj ~ 5300 A stand out in the 2D plot of WRMS 
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Fig. 7. Results from 10-fold cross-validation on maximum-light spectra. From top to bottom: H only; {x\,'R)\ {c,K)\ {x\,c,K). The left column 
is color-coded according to the weighted rms of prediction Hubble residuals, while the right column coiTesponds to the absolute Pearson cross- 
correlation coefficient of the correction terms with uncorrected Hubble residuals. 



residuals, and all the top color-corrected ratios for {x\ , c, in- 
clude a wavelength bin in that region (which corresponds to the 
absorption trough of the S II /15454 line). 



4.3.3. Comparison with Bailey et al. (2009) 

We confirm the basic result of Bailey et al. (2009) using an inde- 
pendent sample and a different cross-validation method: the use 
of a flux ratio alone or in combination with a color parameter re- 
sults in a Hubble diagram with lower scatter when compared to 
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the standard (xi, c) model. Using a flux ratio alone, Bailey et al. 
(2009) find '^(6420/4430) as their most highly-ranked ratio, 
while we find -^(6630/4400) [see Table 2]. The wavelength bins 
are almost identical, and in any case '??(6420/4430) is amongst 
our top 5 ratios. For this ratio we find j - -3.40 + 0. 10, in agree- 
ment with J = -3.5 ± 0.2 found by Bailey et al. (2009)1 

The other four flux ratios given by Bailey et al. (2009) [their 
Table 1] are not part of our top-5 %. For two of these ratios 
the reason is trivial: they include wavelength bins redder than 
7100 A, not covered by most of our spectra. The other two flux 
ratios ['^(6420/4170) and ??(6420/5120)] lead to differences 
< 5% on the Hubble diagram residual scatter with respect to 
the standard (xi, c) model according to Bailey et al. (2009) [cr = 
0. 1 66 + 0.0 1 6 mag for '??(6420 /4 1 70) and cr = 0. 1 54 + 0.0 1 5 mag 
for 9?(6420/5120), cf. 0.161 + 0.015 mag for (xi,c)], and they 
rank 29* and 658* in our study, respectively. This discrepancy 
is in part due to the selection method: Bailey et al. (2009) select 
their best ratios based on cross-correlation coefficients with un- 
corrected magnitudes, while we select them based on the intrin- 
sic prediction error from cross- vahdated Hubble diagram resid- 
uals. However, ranking our ratios using the same method as 
Bailey et al. (2009) does not resolve the discrepancy. It is pos- 
sible that Bailey et al. (2009) are sensitive to their exact choice 
of training and vaUdation samples, where we have randomized 
the approach. We note however that the impact on the weighted 
rms of prediction residuals is statistically indistinguishable for 
many flux ratios given our sample size (e.g. error on WRMS 
~ 0.03 mag cf. differences of < 0.01 mag in WRMS for the top 5 
flux ratios; see Table 2), so that the exact ranking of flux ratios is 
not well determined and subject to revisions from small changes 
in the input data. 

Using both a color-corrected flux ratio and the SALr2 
color parameter decreases the residual scatter further, as 
found by Bailey et al. (2009). Using the set of predictors 
[c,'??''(6420/5290)] leads to ~ 15% lower WRMS with respect 
to (xi,c) [WRMS = 0.175 + 0.025 mag cf. 0.204 + 0.029 mag], 
and to ~ 20% lower o-pjed [o-pred = 0.148 + 0.029 mag cf. 
0.181 ± 0.032 mag] at \Aa significance based on the difference 
in intrinsic prediction error, t^Xi,c- None of the color-corrected 
flux ratios listed by Bailey et al. (2009) [their Table 1] are part 
of our five highest-ranked although our top ratios are formed 
with almost the same wavelength bins [!R'^(6420/5290) in this 
paper; '^'^(6420/5 190) in Bailey et al. (2009)]. The oflier color- 
corrected ratios in Bailey et al. (2009) rank well below in our 
study, whether we select the best 'J^ according to the resulting 
Hubble residual scatter or the cross-correlation of (-jSc + YK') 
with uncorrected residuals. The same caveats apply here as when 
selecting the best uncorrected flux ratios (see previous para- 
graph), although the 1^ measurement is probably even more sen- 
sitive to the relative flux calibration accuracy of the spectra. 

We also cross-checked the results of Bailey et al. (2009) 
by simply validating their best flux ratios on our entire SN la 
sample. The results are displayed in Table 3, where we give the 
weighted rms of Hubble residuals from a simultaneous fit to the 
entire SN la sample (as done by Bailey et al. 2009), as opposed 
to prediction residuals under cross-vaUdation. For all flux ratios 
(both % and 'R') in Table 3, our own best-fit j agrees within the 
Icr errors with that found by Bailey et al. (2009) [noted 7(B09) 
in Table 3] , although we have systematically larger errors. We 
note that most of the top ratios reported by Bailey et al. (2009) 



^ In fact Bailey et al. (2009) find y = +3.5 + 0.2, but this is due to a 
typo in their equation for the distance modulus: y'R really appears as a 
negative term in their paper (S. Bailey 2010, private communication). 



Table 3. Validation of top 5 flux ratios at maximum light from Bailey et al. (2009) [noted 
B09] 



ix -ty r y(B09) WRMS" A^^^^ 

n 

6420 4430 -3.40 ±0.32 -3.5 ± 0.2 0.184 ±0.026 -0.007 ± 0.031 (0.2o-) 
6420 4170 -4.43 ±0.41 -4.9 ± 0.2 0.197 ±0.028 0.008 ± 0.028 (0.3tr) 
7720* 4370 ■■■ 7.3 ± 0.3 

6420 5120 -4.34 ±0.42 -4.7 ± 0.3 0.237 ± 0.032 0.052 ± 0.034 (l.Scr) 

7280* 3980 7.9 ± 0.3 — 

(c,«') 

6420 5190 -2.20 ±0.57 -3.5 ± 0.3 0.176 ± 0.025 -0.012 ± 0.016 (0.8a-) 
5770 6420 0.67 ±0.22 1.4 ±0.1 0.193 ± 0.027 0.005 ± 0.011 (0.5a-) 
6420 5360 -1.71 ±0.43 -2.3 ±0.2 0.169 ± 0.024 -0.019 ± 0.016 (1.2a-) 
6760 6420 1.79 ±1.09 4.2 ±0.5 0.217 ± 0.032 0.030 ± 0.025 (1.2a-) 
6420 4430 -2.56 ±0.63 -3.2 ±0.3 0.168 ± 0.024 -0.027 ± 0.023 (1.2a-) 

(^i.c) 

■■■ ■■■ 0.194 ±0.027 



Notes. Weighted rms of Hubble residuals from a simultaneous fit 
to the entire SN la sample (as done by Bailey et al. 2009), as opposed 
to prediction residuals under cross-validation. As explained in § 2.3, 
the weighted rms of prediction Hubble residuals is a more realistic esti- 
mate of the accuracy of a given model in measuring distances to SN la. 

Wavelength bins redder than 7 100 A, not covered by most of our 
spectra. 

lead to no significant improvement over (xi,c), and even leads 
to slightly worse results for some ratios [e.g. ??(6420/5120) re- 
sults in WRMS = 0.237 + 0.032mag cf. 0.194 + 0.027 mag 
for (xi,c)]. A closer look at Table 1 of Bailey et al. (2009) 
shows that this is also the case in their paper: for the 'R-only 
model, only one ratio out of five, namely ??(6420/4430), re- 
sults in a lower Hubble diagram residual scatter. The other four 
are either consistent with no improvement [??(7720/4370) and 
??(6420/5120)], or yield slightly worse results ['??(6420/4170) 
and ??(7280/3980)]. Again, this results from the way Bailey 
et al. (2009) selected their best ratios, based on the correlation 
with uncorrected Hubble residuals. 

We cannot directly compare the resulting scatter in Hubble 
diagram residuals with those reported in Table 1 of Bailey 
et al. (2009). First, they use the sample standard deviation (cr), 
whereas we use the weighted rms (see § 2.4). Second, the scat- 
ter they find for the standard (xi , c) model is significantly lower 
than ours. We have refit the data presented in Table 1 of Bailey 
et al. (2009) to derive the weighted rms of Hubble residuals 
for the (xi,c) model from their sample, and find WRMS = 
0.148±0.014mag, which is almostO.05 mag smaller when com- 
pared to our sample (0.194 + 0.027 mag). This difference in the 
Hubble residual scatter between the SN-Factory and CfA sam- 
ples is consistent with the difference found amongst other nearby 
SN la samples by Hicken et al. (2009b). 

Interestingly, using the WRMS statistic as opposed to the 
sample standard deviation results in a smaller difference in resid- 
ual scatter between the 'R-only and (xi , c) models. Using our own 
fits of the data presented in Table 1 of Bailey et al. (2009), we 
find WRMS = 0.131 ± 0.014mag for •R(6420/4430), i.e. ~ 11% 
smaller scatter when compared to (xi , c), where the difference 
between the two models is ~ 20% when considering the sample 
standard deviation. 

4.3.4. Comparison with 2D models 

We use synthetic spectra based on a recent 2D survey of delayed- 
detonation SN la models by Kasen et al. (2009) to investigate the 
physical origin of the high correlation between several flux ratios 
and uncorrected SN la magnitudes. These models were found to 
reproduce the empirical relation between peak fi-band magni- 
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tude and post-maximum decline rate. A more detailed compari- 
son of SN la data with these models wiU be presented elsewhere. 

We measured flux ratios in the same manner as we did for our 
data, and computed Pearson correlation coeflicients with (uncor- 
rected) absolute magnitudes synthesized directly from the spec- 
tra. The 2D correlation map is shown in Fig. 9 (left panel), along- 
side the same map derived from the CfA SN la sample (right 
panel). At first glance, the two maps appear similar, with two 
large ~ lOOOA-wide "bands" of flux ratios with strong corre- 
lations with uncorrected magnitudes, for Ax(Ay) ^ 6200 A and 
Ay(Ax) ^ 6000 A, although the correlations are even stronger 
in the models (several flux ratios have absolute Pearson corre- 
lation coefficients |r| > 0.95, where there are none in the data). 
A closer look reveals some important difl'erences, the models 
having strong correlations for 6000 A< Ax(Ay) < 6200 A and 
AyiAx) ^ 6200 A that are not present in the data. The same ap- 
plies to the regions with coordinates A-xiAy) « 4200 A. These 
differences are significant and illustrate the potential for such 
comparisons to impose strong constraints on SN la models. 

In Fig. 10 we show the correlation of uncorrected abso- 
lute rest-frame 5-band magnitudes (Mb) with our highest-ranked 
flux ratio '^(6630/4400), both from flie 2D models and CfA 
data, where we have used the redshift-based distance for the 
latter. The vertical offset is arbitrary and solely depends on 
the normalization adopted for the data, which we have chosen 
for sake of clarity. There are 1320 model points, each corre- 
sponding to one of the 44 2D delayed-detonation models of 
Kasen et al. (2009) viewed from one of 30 different viewing an- 
gles. The linear fits shown in Fig. 10 are done over the range 
0.25 < '^(6630/4400) < 0.50, where the models and data over- 
lap. For the data this is equivalent to excluding the three most 
highly-reddened SN la (open circles), for which the host-galaxy 
visual extinction Ay was determined based on light-curve fits 
wifli MLCS2k2 (Jha et al. 2007). This is justified since no red- 
dening by dust is applied to the models. The slope of the relation 
between Mb and '^(6630/4400) is significantly steeper for the 
models (F = 7.38 ± 0.13) flian for the data (F = 4.76 + 1 .04), and 
the correlation is much stronger (r = 0.89 cf. 0.69 for the data). 
This is not surprising since the data are subject to random mea- 
surement and pecuUar velocity errors, which degrade the corre- 
lation. Including models for which '??(6630/4400) < 0.25 softens 
the slope to F = 5.71 + 0.03 and results in a stronger correlation 
(r = 0.97), while including data with Ay > 0.45 mag results in 
F = 4.43 ± 0.47 and a much stronger correlation r = 0.92. This 
last value for F can be compared with the y fitting parameter for 
this same flux ratio (y = -4.42 + 0.09; see Table 2), although 
the latter is based on a formal cross-vahdation procedure and 
the opposite sign is a consequence of the convention when using 
the flux ratio to predict SN la distances. As noted in § 4.3.2, the 
correlation of Mb with '^(6630/4400) is largely biased by the 
minority of highly-reddened SN la. 

The models yield values for ??(663 0/4400) ranging between 
~ 0.12 and ~ 0.44, all due to intrinsic color variations. Since 
these models reproduce the relation between Mb and post- 
maximum decline rate of Phillips ( 1 993), they confirm the intrin- 
sic nature of the correlation between '??(6630/4400) and {xi,c} 
shown in Fig. 6 (upper left panel). 

The wavelength bins Ax - 6630 A and Ay = 4400 A are 
close to the central wavelengths of the standard R and B broad- 
band filters, hence ^(6630/4400) is a rough measure of the B-R 
color at B-band maximum. The 2D models of Fig. 10 indicate 
that a large part of the variation in '??(6630/4400) seen in the data 
is due to intrinsic variations in B-R color. Reddening in the host 
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Fig. 10. Absolute rest-frame 5-band magnitude (Mb) vs. flux ratio 
'/?(6630/4400) at maximum light in 2D SN la models of Kasen et al. 
(2009) [.small dots], and in data from the CfA SN la sample (Ay < 
0.45 mag: filled circle.^. Ay > 0.45 mag: open circles). The dotted 
and dashed lines are linear fits to the models and data, respectively, 
where models for which ??(6630/4400) > 0.25 and data for which 
Av > 0.45 mag have been excluded from the fit. The slope (F) and 
Pearson correlation coefficient (r) are indicated for both cases. The 
data have been offset vertically for clarity. Including models for which 
;R(6630/4400) < 0.25 results in F = 5.71 ± 0.03 and r = 0.97, while in- 
cluding data with Av > 0.45 mag results in F = 4.43±0.47 and r = 0.92. 
The arrows indicate approximate reddening vectors for different values 
ofRy. 

galaxy is then needed to explain values of !R(6630/4400) > 0.4, 
while at lower values it is challenging at best to discriminate be- 
tween the effects of intrinsic color variations and extinction by 
dust, since both affect '^(6630/4400) in the same manner, as il- 
lustrated by the reddening vectors in Fig. 10 [they are really red- 
dening curves, cf. Fig. 3, but the behavior is almost linear over 
this smafl range in '^(6630/4400)]. 

The models also give a physical explanation for the corre- 
lation of :R(6630/4400) with absolute magnitude. Indeed, the 
variation of this ratio is largely caused by spectroscopic varia- 
tions around 4400 A, a region dominated by lines of Fe 11 and 
Fe 111, with contributions from Mg 11 (Ti 11 provides an impor- 
tant source of opacity for the least luminous SN la), while the 
region around 6630 A has little intrinsic variation (this was noted 
by Bailey et al. 2009). This translates to a standard deviation of 
peak B-band magnitudes (cr » 0.40 mag) that is almost twice as 
large as the /?-band magnitude (at B maximum; cr » 0.26 mag) 
in the models. The relative contribution of Fe II and Fe III lines 
is related to the temperature of the line-forming regions in the 
SN la ejecta, itself a function of peak luminosity (dimmer SN la 
are generally cooler; see, e.g., Kasen & Woosley 2007). One thus 
expects a large luminosity-dependent spectroscopic variation in 
this wavelength region, although its exact shape and relation to 
temperature remains largely unknown. 

While these models provide useful insights into the physical 
origin of these correlations, a direct comparison with the data 
reveals some of their shortcomings. In Fig. 10 we see that some 
models predict values of the flux ratio '??(6630/4400) < 0.25 
for the most luminous SN la, where the data are limited to val- 
ues greater than this. Our sample includes several SN la at the 
high luminosity end that show no sign of extinction in their 
host galaxies (Ay < 0.05 mag based on light-curve fits with 
MLCS2k2), so the differences are real and point to discrepan- 
cies between the data and the models, some of the latter having 
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Ax [A] 

Fig. 9. Absolute Pearson correlation coefficients of flux ratios at m 
detonation SN la models of Kasen et al. (2009) (left), and in data from 




num light with uncorrected absolute magnitudes Mg in 2D delayed- 
CfA SN la sample (right). 



bluer B -7? colors at Z?-band maximum. This is not surprising, as 
the models explore a larger range of parameter space than is re- 
alized in nature. Comparisons of this sort can then help constrain 
the range of model input parameters. A more detailed compari- 
son of SN la data from the CfA SN program with these models 
will be presented elsewhere. 

4.4. Results on spectra at other ages 

Bailey et al. (2009) restricted their analysis to spectra within 
Af = 2.5 d from B-band maximum. In this section we consider 
flux ratios measured on spectra at other ages. We impose the 
same cuts on relative flux calibration accuracy (\A(B - V)\ < 
0.1 mag), SALT2 color (c < 0.5), redshift (z > 0.015), and age 
range (At = 2.5 d) as those used for the maximum-light spectra 
in the previous section. We consider all ages between t = -2.5 d 
and t = +7.5 d, in steps of 2.5 d (for ages earlier than -2.5 d or 
later than -1-7.5 d the number of SN la with spectra that satisfy 
our cuts falls below 20, and we do not trust the results). We re- 
port the best ratio at each age in Table 4, for both the '7?-only and 
(c, •/?') models. 

The 2D maps of Hubble diagram residual scatter for spec- 
tra at ages t = -2.5, -t-0, +5, and -1-7.5 d are shown in Fig. 11. 
As was the case for maximum-light spectra, adding the SALT2 
xi parameter leads to slightly degraded results when compared 
with 'R alone, so we do not show plots for (xi,'R) in Fig. 11. 
Moreover, we do not show results for (xi , c, 'R'^) since the best 
color-corrected flux ratios in this case do not result in a sig- 
nificant improvement over (c, 'R'^). At all the ages we consider 
here, the set of predictors (c, 'R'^) results in lower weighted rms 
of prediction residuals than (xi,c), although the significance 
of the difference is < 2cr and is lower for f > -1-5 d than for 
-2.5 < f < -1-2.5 d. A flux ratio by itself only leads to an im- 
provement over (xi,c) near maximum light (-2.5 < f < -1-2.5 d). 
As was the case at maximum light, there is a positive intrinsic 
correlation in prediction error between all the distance prediction 
models that include a flux ratio and (xi,c) [0.4 < pi,,c < 0.8]. 

The best set of predictors overall in the age range -2.5 < f < 
-1-7.5 d is [c, :R''(4610/4260)] at t = -2.5 d. The weighted rms of 
prediction residuals is reduced by ~ 30% with respect to (xi , c) 
[WRMS = 0.143 +0.020 mag], and the intrinsic prediction error 
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Fig. 12. Correlation between (46 10/4260) at t = -2.5 d and the 
SALT2 fit parameters (xi,c). The Pearson coefficient of the correlation 
with color drops to r = 0.03 if we ignore the two points at c > 0.4. 



by ~ 40% (cTpi-ed = 0.106 + 0.028 mag), the significance of the 
difference being ~ 2cr (A^^, ^ = -0.081 + 0.037 mag). We show 
the correlation of (46 10/4260) at t = -2.5 d and the SALT2 
parameters (xi,c) in Fig. 12. This color-corrected ratio is mildly 
correlated with xi [r - 0.61; this drops slightly to r = 0.47 if 
we ignore the two points at '^'^(4610/4260) > 1.3] and uncorre- 
cted with color (r - -0.40; this drops to r = 0.03 if we ignore 
the two points at c > 0.4). Interestingly, the wavelength bins 
that constitute this ratio are part of the two prominent spectral 
absorption features, predominantly due to iron-group elements, 
that were found to vary intrinsically between SN la based on the 
2D models discussed in § 4.3.4. 



4.5. Results using two flux ratios 

In this section we consider corrections using a linear combina- 
tion of two flux ratios as follows: 



ju = niB - M + jiRi + J2R2 



(14) 
(15) 



where the latter equation includes an additional correction due to 
color, hence the use of color-corrected ratios For both cases. 



13 



S. Blondin et al.: Do spectra improve distance measurements of Type la supemovae? 



Table 4. Top flux ratio at ages -2.5 < t < +7.5 d from 10-fold CV 



Rank 


Ax 


Ay 


r 


WRMS 


Cpred 








K 


-2.5 


6540 


4580 


-6.09 + 0.11 


0.182 


+ 0.025 


0.151 +0.031 


0.70 + 0.13 


-0.032 + 0.028 (l.lcr) 


24 


+0.0 


6630 


4400 


-4.37 ± 0.09 


0.189 


+ 0.026 


0.163 +0.030 


0.80 + 0.09 


-0.018 + 0.025 (0.7cr) 


26 


+2.5 


6630 


4040 


-3.51+0.09 


0.203 


+ 0.027 


0.171 +0.033 


0.63 + 0.13 


-0.012 + 0.033 (0.4cr) 


26 


+5.0 


6590 


4490 


-4.69 + 0.12 


0.225 


+ 0.031 


0.203 + 0.034 


0.36 + 0.20 


0.022 + 0.041 (0.5cr) 


26 


+7.5 


6590 


4890 


-3.40 ± 0.20 


0.251 


+ 0.035 


0.229 + 0.039 


0.47 + 0.18 


0.044 + 0.043(1.00") 


25 




-2.5 


4610 


4260 


2.19 + 0.14 


0.143 


+ 0.020 


0.106 + 0.028 


0.43 + 0.23 


-0.081 +0.037 (2.20-) 


24 


+0.0 


6420 


5290 


-1.75+0.10 


0.175 


+ 0.025 


0.148 + 0.029 


0.80 + 0.09 


-0.032 + 0.023 (1.4cr) 


26 


+2.5 


5550 


6630 


1.09 + 0.08 


0.169 


+ 0.022 


0.133 +0.030 


0.79 + 0.09 


-0.049 + 0.027 (1.8(j) 


26 


+5.0 


6540 


5580 


-5.18 + 0.48 


0.194 


+ 0.026 


0.166 + 0.031 


0.57 + 0.16 


-0.014 + 0.034 (0.4O-) 


26 


+7.5 


6460 


5510 


-1.85 + 0.13 


0.200 ± 0.028 


0.173 + 0.033 


0.79 ± 0.09 


-0.007 + 0.027 (0.3(r) 


25 
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Fig. 13. Correlation between tlie highest-ranked (^2, t^) at maximum 
light and the SALT2 fit parameters {jci,c), and the highest-ranked 



we fix f6^^ to the highest-ranked single flux ratio [e.g., at max- 
imum light: (6630/4400) and (6420/5290); see Table 2], 
but leave both (71,72) as free parameters (i.e. we do not set 71 
equal to 7 found in the single flux ratio case). 

The results for the top five secondary flux ratios at max- 
imum light are displayed in Table 5. In all cases, includ- 
ing a second flux ratio further reduces the weighted rms 
of prediction residuals by * 15-20% [WRMS = 0.162 ± 
0.022 mag for ^2(5160/5290); WRMS = 0.151 ± 0.021 mag 
for '^2(5690/5550)] with respect to the single flux ratio case 
[WRMS = 0.189 + 0.026 mag for ^(6630/4400); WRMS = 
0.175 ± 0.025 mag for (6420/5290)]. Again flie significance 
of the improvement (< 2cr) is difficult to gauge given our 
sample size, this despite the fact that 72 is significantly dif- 
ferent from zero in all cases. Our best secondary flux ratios, 
:R2(5 160/5290) and ^^^(5690/5550), are uncorrelated with the 
SALT2 fit parameters {x\,c) and with the highest-ranked pri- 
mary ratios (6630/4400) and (6420/5290) [see Fig. 13], 
and hence provide independent information on SN la luminos- 
ity. 

The secondary flux ratios 'R2 listed in Table 5 have a 
much smaller wavelength baseline than the primary ratios Hi 

and "^^j (see Table 2). The highest-ranked secondary ratios 
have a 130A and 140A baseline, respectively. These ratios 



do not measure SN la colors: they measure small-scale intrin- 
sic spectroscopic variations. Interestingly, all the wavelength 
bins that form these secondary ratios are clustered around the 
S II /1/15454,5640 doublet and the iron-group-dominated absorp- 
tion complex Fe 11/14800 mentioned in § 4.4 (see also § 4.3.4). 
As was the case for a single flux ratio, the results using a sec- 
ondary flux ratio are similar for -2.5 < t < +2.5 d, and tend to 
be worse at later ages. We choose not to discuss them further. 

In a recent paper"^, Yu et al. (2009) have searched for flux ra- 
tio pairs that minimize Hubble diagram residuals (with no color 
correction), and find several such pairs which achieve a stan- 
dard deviation cr < O.lOmag at ages between -3d and +12d 
from maximum light. We have validated the flux ratio pairs re- 
ported in their Table 4 and find that none of them lead to an 
improvement compared with the standard (xi,c) model. There 
could be several reasons for this disagreement: Yu et al. (2009) 
do not use a cross-validation procedure and their sample size 
(anywhere between 17 and 24 SN la depending on the age and 
flux ratio pair considered) suggests they may be overfitting a 
small sample. Moreover, they use color-corrected flux ratios (ac- 
tually corrected for the host-galaxy extinction, Ay, as opposed to 
the SALT2 color parameter) but do not include a color parameter 
in their equation to correct for the SN la magnitudes. When we 
use the S ALT2 color parameter in addition to the same flux ratio 
pairs as Yu et al. (2009), several pairs indeed lead to an improve- 
ment over the standard {x\,c) model, but not over the correction 
using a single flux ratio (or a single flux ratio in combination 
with color) we find in this paper. Last, Yu et al. (2009) do not 
impose a redshift cut on their sample: only 14 out of 38 SN la 
are at redshifts z > 0.015, and 5 are at redshifts z < 0.005, where 
the magnitude error due to peculiar velocities is CTptc > 0.4 mag. 
It is unclear why this should lead to a lower scatter in Hubble 
residuals (on the contrary one expects an increased scatter), but 
it certainly impacts their analysis. 



5. Other spectroscopic indicators 

In this section we consider other spectroscopic indicators, 
mostly related to spectral line profile morphology. Some of these 
indicators are also flux ratios, but the wavelengths correspond to 
precise locations of absorption troughs or emission peaks in the 
SN la spectrum, as opposed to the "blind" approach of comput- 
ing flux ratios from all possible wavelength bins with no a priori 
physical motivation. We use the same approach as for the flux 



At the time of writing, the paper by Yu et al. (2009) has not been 
accepted for publication. Here we refer to the 2'"' version of their paper, 
dated 30"" June 2010. 
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Fig. 11. Results from 10-fold cross-validation on spectra at / = -2.5, +0, +5, +7.5 d. (from top to bottom), color-coded according to the weighted 
rms of prediction Hubble residuals. The left column is corresponds to H only, while the right column corresponds to the (c, 'R'^) model. 



ratios, i.e. we consider models for predicting SN la distances 
which include a spectroscopic indicator, possibly in combina- 
tion with a light-curve parameter (cf. Eq. 1), and we validate 
each model using /T-fold cross-validation (and present results 
for = 10 in this section). 



5.1. Measurements 

We divide the SN la spectrum into several "features", each la- 
beled according to the strongest line in that wavelength range. 
Figure 14 shows the seven features we consider in this pa- 
per, from CaIU3945 in the blue to SiIU6355 in the red. The 
wavelengths associated with each feature correspond to the gf- 
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Table 5. Top 5 secondary flux ratios at maximum light from 10-fold CV on 26 SN la 



Rank Ax 


Ay 


71 




72 


WRMS 




P-v, 


.C 






[7?! (6630/4400), 


1 5160 


5290 


-5.05 ± 


0.17 


-2.24 + 0.22 


0.162 


+ 0.022 


0.127 


+ 0.028 


0.32 + 


0.23 


-0.052 


+ 0.038 (1.40-) 


2 5290 


5160 


-5.08 ± 


0.18 


4.78+0.50 


0.167 


+ 0.022 


0.133 


+ 0.029 


0.31 + 


0.24 


-0.046 


+ 0.036 (1.30-) 


3 5120 


5290 


-5.19 ± 


0.11 


-1.75+0.20 


0.172 


+ 0.024 


0.141 


+ 0.029 


0.41 + 


0.21 


-0.038 


+ 0.036 (1.10") 


4 5690 


5360 


-3.93 ± 


0.12 


-2.06 + 0.21 


0.172 


+ 0.024 


0.141 


+ 0.029 


0.46 + 


0.20 


-0.038 


+ 0.032 (1.20") 


5 5690 


5290 


-4.14 ± 


0.11 


-1.26 + 0.13 


0.172 


+ 0.024 


0.142 


+ 0.029 


0.48 + 0.19 


-0.038 


+ 0.032 (1.20") 


[C;?'; (6420/5290), 7?:;] 


1 5690 


5550 


-1.69 + 0.09 


-1.93+0.21 


0.151 


+ 0.021 


0.115 


+ 0.027 


0.58 + 


0.17 


-0.063 


+ 0.030 (2.I0-) 


2 5690 


5960 


-2.34 ± 


0.11 


-1.90 + 0.23 


0.152 


+ 0.020 


0.116 


+ 0.028 


0.42 + 


0.21 


-0.064 


+ 0.035 (I.80-) 


3 5960 


5690 


-2.35 ± 


0.11 


2.05 + 0.24 


0.153 


+ 0.021 


0.117 


+ 0.028 


0.43 + 


0.21 


-0.063 


+ 0.036 (1.80") 


4 5550 


5690 


-1.70 ± 


0.09 


1.72 + 0.20 


0.152 


+ 0.021 


0.118 


+ 0.027 


0.58 + 


0.17 


-0.061 


+ 0.031 (2.O0") 


5 4890 


5690 


-1.25 ± 


0.10 


1.22 + 0.24 


0.158 


+ 0.021 


0.122 


+ 0.028 


0.62 + 


0.15 


-0.052 


+ 0.026 (2.O0-) 
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Fig. 14. Wavelength bounds of spectroscopic features for which we 
measured the various indicators shown in Fig. 15, illustrated using the 
maximum-hght spectrum of SN 2006ax. 



weighted mean wavelength of the different atomic transitions for 
the specified ion (e.g. 3945 A for the Call H & K lines), except 
for the two large features dominated by Fe II lines, where the 
wavelengths denote the approximate location of the deepest ab- 
sorption (~ 4300 A and ~ 4800 A). The FeIU4300 feature also 
includes contributions from Mgll, and possibly Felll for the 
most luminous SN la. For the faintest, 1991bg-like SN la. Till 
constitutes a dominating source of opacity in this wavelength re- 
gion. Since we do not include 1991bg-like SN la in our analysis, 
however, we do not present measurements for Ti II. 

The various spectroscopic indicators we consider are illus- 
trated in Fig. 15, based on the Si 11/16355 line profile in the 
spectrum of SN 2005ki at t - +ld from maximum light. We 
first smooth the spectrum using the inverse-variance Gaussian 
filter of Blondin et al. (2006) with a smoothing factor 0.001 < 
dA/A < 0.01 determined based on ax^ test using flux errors from 
the variance spectra (Fig. 15; thick line). The smoothed spectrum 
makes it easier to define wavelength locations of local flux max- 
ima on either side of the absorption component of the P Cygni 
profile (/Ibiue and /Ipeak), as well as the location of maximum ab- 
sorption (/iabs)- The wavelengths /labs and /Ipeak are then used to 
define the absorption and peak velocities, respectively (vabs and 
Vpeak), using the relativistic Doppler formula (see also Blondin 
et al. 2006). We also measure the heights of the local maxi- 
mum (/ibiue and /ipeak), and define a pseudo-continuum between 



i 

Pi 




Fig. 15. Definition of the main spectroscopic indicators used in this pa- 
per, here illustrated using the Si 11 /i6355 line profile in the spectrum of 
SN 2005ki at ? = -1- 1 d. The right panel shows the pseudo-continuum 
(dashed line), as well as the wavelength locations of the blue and red 
emission peaks (/Ibiue and /ipc^k) and their respective heights (/ibiue and 
/ipcak). The wavelength of maximum absorption (/labs) serves to define 
the absorption velocity, v^bs. The peak velocity Vpeak is defined analo- 
gously. The left panel shows the same line profile normalized to the 
pseudo-continuum, and serves to define the (relative) absorption depth 
(rfiibs)j FWHM, and pseudo-equivalent width (pEW; shaded gray re- 
gion). In both panels, the thick line corresponds to the smoothed flux, 
where we have used the inverse-variance weighted Gaussian filter of 
Blondin et al. (2006) with a smoothing factor dA/A = 0.005. 



them. These latter quantities are measured on the original, un- 
smoothed spectrum. Division by this pseudo-continuum enables 
us to measure the relative absorption depth (c/abs) and full-width 
at half-maximum (FWHM) of the absorption component, as well 
as its pseudo-equivalent width (pEW; defined analogously to the 
equivalent width used by stellar spectroscopists for abundance 
determinations, but without the physical basis, hence "pseudo" 
EW; Fig. 15; right panel). 

We measure these quantities for all the features presented in 
Fig. 14, except for the complex Fe II /14300 and Fe II /14800 fea- 
tures for which we only consider the pseudo-equivalent width. 
The error on each measured quantity includes errors due to red- 
shift, relative flux calibration, host-galaxy extinction and con- 
tamination, and of course the flux error. We only consider mea- 
surements for which the mean S/N over the entire feature is 
greater than 5 per A, and require a minimum of 20 SN la with 
valid measurements. Note that we do not impose cuts on relative 
flux calibration accuracy or SALT2 color, as was the case for the 
flux ratios, since these quantities are mostly local measurements 
which are far less sensitive to the overall SED. 
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We also consider various spectroscopic ratios, which were 
found to correlate with absolute magnitude, defined below: 



«(Ca) = 



^(CaS) = 



^(Si) 



KiSiS) 



??(SiSS) 



??(S,Si) = 



??(Si,Fe) 



Vak(CalU3945) 
/!biue(CaIU3945) 

J3887 ^^^'^ 

t/ab^(Si 11 .15972) 
(/abs(Si 11 .16355) 
/!peak(Si 1U6355) 
/!pe„k(S 1U5640) 

^5700 
J55OO 



FidA 



J6200 

pEW(S 11/1^5454,5640) 

pEW(Si 11/15972) 
pEW(Si 11/15972) 
pEW(Fe 11/14800)' 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



The ratios '??(Ca) and 'RiSi) were both defined by Nugent et al. 
(1995), and found to correlate well with the luminosity decline 
rate parameter AOTi5(fi). To increase the S/N of the HiCa) mea- 
surement, Bongard et al. (2006) introduced the corresponding 
integral flux ratio !R(CaS), also found to correlate with absolute 
magnitude. Using a grid of LTE synthetic spectra to investigate 
the 'R(Si) wavelength region, Bongard et al. (2006) also defined 
a ratio of the red local maximum of Sill /16355 to the red local 
maximum of S II /15640, noted !??(SiS). The corresponding inte- 
gral flux ratio is '^(SiSS), again introduced by Bongard et al. 
(2006) to increase the S/N of the 'R(SiS) measurement. Last, 
Hachinger et al. (2006) measured the absorption velocities and 
pseudo-EW in 28 SN la spectra and found two additional pEW 
ratios, 'R(S, Si) and '^(Si, Fe), that are good indicators of lumi- 
nosity. Note that !??(Ca) and '??(SiS) are in fact flux ratios similar 
to fliose defined by Bailey et al. (2009). 

5.2. Results 

We present our results using the absorption velocity (v^bs; units 
of 10"* km s"'), the fuU-width at half-maximum (FWHM; units 
of 10^ A), the relative absorption depth (Jabs), the pseudo- 
equivalent width (pEW; units of 10^ A), and the various spectro- 
scopic ratios KiX) [Eqs. 16-22] in Tables C. 1-C.5 (Appendix C). 
We do not present results for the peak velocity (Vpeak) as they are 
far worse than for the other indicators. There were not enough 
valid measurements for Ca 11/13945, hence the absence of this 
line in Tables C.1-C.4. We only report results for the bluer ab- 
sorption of the SII doublet (SIU5454) in Tables C.1-C.3, but 
the pseudo-equivalent width is that of the entire doublet (see 
Table C.4). 

Based on the difference in intrinsic prediction error with re- 
spect to the standard model which uses the SALT2 fit parameters 
(xi, c), again noted Aj^^^c^ we see from Tables C.1-C.5 that none 
of these spectroscopic indicators alone leads to a lower weighted 
rms of prediction residuals (i.e. A^, _c > 0). At best they are con- 
sistent with no improvement at all [e.g. pEW(Si 11/14 130), for 
which A;c,_c = 0.041 + 0.033 mag)]. The same is true at ages 
other than maximum Ught. 

Nonetheless, several such indicators compete well with 
(xi,c), even leading to small improvements (albeit statisti- 
cally insignificant), but only when combined with SALT2 
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Fig. 16. Correlation between pEW(SiIU4130) and HiSi) at maximum 
light and the SALT2 fit parameters (xi,c), and color-corrected Hubble 
residual. The open circle in the lower panels corresponds to SN 2000dk. 



color [pEW(SiIU4130) and •R(Si)] or in addition to (xi,c) 
[Vabs (Si 11/16355 and c/abs(S 11/15454)]. We discuss these indica- 
tors in the two following sections. 



5.2.1 . Spectroscopic indicators in combination witli SALT2 

color 

When used in combination with the SALT2 color parameter, 
both the pseudo-equivalent widfli of Si 11/14130 and flie 'R(Si) 
spectroscopic ratio compete well with the standard (xi , c) predic- 
tors (A;t,,c = 0.006±0.014magand A;c;,c = -0.007+0.030mag). 
Both indicators are strongly anti-correlated with x\ and uncorre- 
cted with SALT2 color (see Fig. 16, left and middle panels), 
while the correlation with color-corrected Hubble residual is 
more pronounced for «(Si) [r = 0.63] than for pEW(Sin/14 130) 
\r = 0.35]. In a sense, both indicators act to replace the 
light-curve width parameter xi. The anti-correlation of KiSi) 
with xi has been recovered by several authors since its pub- 
lication by Nugent et al. (1995), while the relation between 
pEW(SiII /14130) and Ught-curve shape has more recently been 
mentioned by Arsenijevic et al. (2008) and Walker et al. (2010). 

We show the Hubble residuals obtained when using 
pEW(SiIU4130) and ??(Si) in combination with SALT2 color 
in Fig. 17, where we also show the residuals from the stan- 
dard (xi,c) model. The subluminous (but not 1991bg-like) 
SN 2000dk stands out as a < 2cr outher for [c,'R(Si)], while 
this is not the case for (xi,c) [the point corresponding to 
SN 2000dk is highlighted in bofli Figs. 16 and 17]. This 
single SN contributes a large fraction of the residual scat- 
ter [WRMS (incl.OOdk) = 0.190 + 0.025 mag cf. 0.196 + 
0.027 mag for (xi,c)], and excluding it from the sample leads 
to a ~ 10% decrease in the weighted rms of prediction Hubble 
residuals, resulting in a ~ 15% improvement over (xi,c) 
[WRMS (excl.OOdk) = 0.171 ±0.028 mag cf. 0.197 ±0.028 mag 
for (xi,c)]. 

5.2.2. Spectroscopic indicators in addition to the SALT2 fit 
parameters (xi,c) 

When used in addition to the standard SALT2 fit parameters 
(xi,c), both the absorption velocity of Si 11/16355 and the rel- 
ative absorption depth of S 11/15454 result in a < 10% de- 
crease in the weighted rms of prediction residuals (Aj, ^ = 



-0.020 ±0.019 mag and 



-0.022 ± 0.030 mag), although 
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Fig. 18. Correlation between Vabs(Sin/l6355) and c/abs(S 11/15454) at 
maximum light and the SALT2 fit parameters {xi,c), and {xi,c)- 
corrected Hubble residual. 
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Fig. 17. Hubble diagram residuals for pEW(Si H ^4130) [top] and ??(Si) 
[bottom] at maximum light. In each case we show the Hubble residu- 
als obtained using SALT2 color and the spectroscopic indicator (up- 
per panels), and using the standard SALT2 fit parameters (xi,c) (lower 
panels). We also indicate the weighted rms of Hubble residuals (gray 
highlighted region). For the !R(Si) spectroscopic indicator, we report the 
weighted rms both including and excluding SN 2000dk (open circle). 



the apparent improvement for (iabs(S H /15454) is due to the fact 
that the weighted rms of prediction residuals for the (xi,c) 
model is somewhat larger for this particular sample (WRMS 
(1.221 ± 0.031 mag). Both these spectroscopic indicators are un- 
correlated with xi or color [Fig. 18; the apparent correlation of 
^abs(Si 11/16355) with color (r - 0.59) is destroyed if we ignore 
the one point at c » 0.2], and thus provide additional information 
independent of light-curve shape or color. The correlation with 
(x] , c)-corrected Hubble residuals (Fig. 18, right panels) is only 
modest (|r| 0.40 for both indicators), and should be reviewed 
as more data become publicly available. 

We show the Hubble residuals obtained when using 
Vabs(Si n /16355) and Jabs(S 11 A5454) in addition to the standard 
(xi,c) predictors in Fig. 19. One clearly sees from these dia- 
grams that the impact of the additional spectroscopic indicator 
is fairly small, as the sign and magnitude of the residuals are al- 
most the same for (xi , c, S) and (xi , c). This is further confirmed 
by looking up the value for the intrinsic correlation in prediction 
error for both indicators in Tables C.l andC.3:pf| c = 0.83+0.06 
for Vabs(Si II /i6355) and p^,,c = 0.75 ±0.10 for d^US II /15454). 



5.2.3. Results using multiple indicators 

We have also considered models involving a linear combination 
of two spectroscopic indicators (i.e. // = niB-M+jiSi +y2S2) or 
a ratio of two indicators (i.e. = ms-M+ySi /S2), also for cases 
including the SALT2 fit parameters (xi,c). No combination of 
two of these spectroscopic indicators leads to an improvement 
over the single indicator case, regardless of the age considered. 



6. Discussion: do SN la spectra really help? 

The central question this paper addresses is whether spectra 
yield useful information to predict distances to SN la better 
than light-curve width and color alone. The answer to this ques- 
tion can have a significant impact on the way future SN la 
surveys are planned, namely whether or not they should in- 
clude spectroscopic (or spectro-photometric) capabilities. This 
has been (and remains!) an active area of discussion for propos- 
als for space-borne missions within the framework of the Dark 
Energy Task Force (Albrecht et al. 2009) or the US Astronomy 
& Astrophysics Decadal Survey^. 

Of all the spectroscopic indicators considered in this paper, 
the concept of flux ratio introduced by Bailey et al. (2009) ap- 
pears to be the most promising, yielding up to ~ 30% lower 
Hubble residual scatter than when using the standard light-curve 
parameters. However, given the limited sizes of the SN-Factory 
(58 SN la) and CfA (26 SN la) samples on which the method 
has been appUed, the results are at best statistically significant 
at the < 2cr level, and the method should be validated on much 
larger samples. It should be noted that the measurement of flux 
ratios requires accurate relative flux calibration, as well as min- 
imal contamination by host-galaxy light. Both requirements im- 
pose strong conditions on future SN la surveys that plan to use 
this method. 

The other spectroscopic indicators we consider in this pa- 
per are intimately linked to fine-profile shapes of specific SN la 
spectral features. One would have hoped that such a physically- 
motivated approach would yield interesting results, but this is 
not the case. At best, these indicators yield < Icr lower resid- 
ual scatter compared with the standard light-curve parameters. 
This is rather disappointing, but also points to potential problems 
with the measurement method we use. It is largely automated, 
but requires some human interaction to ensure the correct local 
maxima used to define the wavelength bounds of each feature 
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Fig. 19. Hubble diagram residuals for Vabs(Si 11/16355) [top], and 
rfabs(S 11/15454) [hotrom] at maximum light. In each case we show the 
Hubble residuals obtained using the spectroscopic indicator in addition 
to the S ALT2 fit parameters (xi , c) {upper panels), and using (xi,c) only 
(lower panels). We also indicate the weighted rms of Hubble residuals 
(gray highlighted region). 



are selected. Moreover, while some indicators (such as the ab- 
sorption velocity Vabs) are largely insensitive to host-galaxy red- 
dening, others (such as the pseudo-EW) are strongly affected. 
Recent unbiased techniques based on wavelet transforms have 
been proposed that are largely insensitive to these measurement 
issues (Wagers et al. 2010), and the present analysis could be 
repeated with such techniques. 

Last, given our spectroscopic data we have focused exclu- 
sively on the optical region, but there appears to be spectroscopic 
indicators that correlate with luminosity in other wavelength re- 
gions (UV: Foley et al. 2008; NIR: Marion et al., in preparation). 
An increased spectroscopic sample at these wavelengths might 
reveal spectroscopic quantities that lead to even more precise 
distances to SN la than optical flux ratios. 

7. Conclusions 

We have investigated the use of spectroscopic indicators which, 
when used alone or in conjunction with light-curve parameters 
(width and color), predict distances to SN la better than when 
using the standard combination of light-curve width and color. 



We have carried our a /T-fold cross-validation analysis on a 
large spectroscopic data set obtained through the Cf A Supernova 
Program. We constructed and implemented maximum likeli- 
hood estimators for the rms intrinsic prediction error of a given 
method, and the intrinsic covariance of prediction errors of dif- 
ferent methods. We used these estimates to compare predictive 
models for SN la distances in a quantitative manner. 

We first considered the spectroscopic flux ratios of Bailey 
et al. (2009), highlighting the importance of an accurate rela- 
tive flux calibration and the failure of this method for highly- 
reddened objects (SALT2 color c > 0.5). At maximum Ught, our 
best single flux ratio !??(6630/4400) from 26 SN la at z > 0.015 
leads to a ~ 10% lower weighted rms of cross- validated predic- 
tion Hubble residuals (WRMS = 0.189 + 0.026 mag) than when 
using the standard SALT2 hght-curve width (xi) and color (c) 
parameters (WRMS = 0.204 ± 0.029 mag), at O.Vcr significance. 
When used in combination with the S ALT2 color parameter, our 
best color-corrected flux ratio '??''(6420/5290) leads to ~ 15% 
lower weighted rms (WRMS - 0.175 + 0.025 mag), at lAcr sig- 
nificance. We thus confirm the use of flux ratios in improving 
distance measurements of SN la magnitudes, although the sig- 
nificance of the difference with respect to the standard purely 
photometric approach is difficult to gauge given our sample size. 
We also point to differences between the best ratios found in 
this paper and those reported by Bailey et al. (2009), in part due 
to the way these ratios are selected: Bailey et al. (2009) select 
their best ratios based on cross-correlation coefficients with un- 
corrected magnitudes, while we directly select them using the 
rms intrinsic error of cross- vaUdated distance predictions in the 
Hubble diagram. 

Comparison of our results with synthetic spectra from a 2D 
survey of delayed-detonation explosion models of Kasen et al. 
(2009) shows that a large part of the variation in our best single 
flux ratio ??(6630/4400) is intrinsic and not due to reddening by 
dust. The correlation of this ratio with SN la magnitudes is due 
to the luminosity-dependent spectroscopic variation in the iron- 
group dominated absorption features around ~ 4300 A. While 
the models confirm the presence of many flux ratios that cor- 
relate strongly with absolute magnitude, significant deviations 
exist with respect to the data. Such deviations can in principle be 
exploited to impose strong constraints on SN la models. 

We extended the analysis of flux ratios to SN la spectra 
at other ages (-2.5 < t < +7.5 d from maximum light (see 
Table 4). The best set of predictors overall in this age range 
is the color-corrected 'R^ (46 10/4260) at f = -2.5 d combined 
with SALT2 color, which leads to ~ 30% lower weighted rms 
of prediction Hubble residuals with respect to (xi, c) [WRMS = 
0.143 ± 0.020 mag], and to ~ 40% lower intrinsic prediction er- 
ror (cTpred = 0.1 06 + 0.028 mag), at ~ 2cr significance. The wave- 
length bins that constitute this ratio are part of the two proininent 
spectral absorption features predominantly due to iron-group el- 
ements, labeled Fe 11/14300 and Fe 11/14800, and which were 
found to vary intrinsically between SN la based on 2D models. 
Flux ratios at f > +5 d fare worse than at maximum light. 

We also considered distance predictions based on two flux 
ratios. We find that the improvement over the standard (xi,c) 
model is at the < 2cr level at best, and tends to be worse for 
ages t > +5 d. At maximum light, our best secondary ratios 
are ??2(5 160/5290) and •^^(5690/5550), whose wavelength bins 
are clustered around the S 11 /1/15454,5640 doublet and the iron- 
group-dominated absorption complex Fe 11/14800. Both ratios 
measure intrinsic small-scale differences between SN la that are 
uncorrelated with light-curve shape or color, and thus provide 
independent information on their luminosity. 
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We also considered spectroscopic indicators associated with 
spectral line-profile morphology: the absorption (and peak) ve- 
locity, the full-width at half maximum, the relative absorption 
depth, the pseudo-equivalent width, as well as other spectro- 
scopic ratios. None of these spectroscopic indicators alone leads 
to a lower weighted rms of prediction Hubble residuals. Only 
when they are combined with SALT2 color do several indicators 
compete well with the standard predictors. Such is the case of the 
Si 11/14130 pseudo-EW and spectroscopic ratio '??(Si). Both in- 
dicators are correlated with xi and act as a replacement to light- 
curve shape in the distance prediction. When used in addition 
to ixi,c), the Si 11/16355 absorption velocity and S 11/15454 rel- 
ative absorption depth lead to a small improvement, albeit sta- 
tistically insignificant. Using a linear combination of two such 
spectroscopic indicators and ratios thereof leads to no further 
improvement, whether at maximum light or at other ages. 

Do spectra improve distance measurements of SN la? Yes, 
but not as much as we had hoped. The statistical framework 
developed here should be applied to an independent and larger 
sample to find out whether the effort of obtaining spectra for 
a cosmological sample will be repaid with better knowledge of 
dark energy. 
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Note added in proof. While this paper was in the final stages of the refereeing 
process, Foley & Kasen (2010) posted a preprint on the arXiv server, in which 
they re-analyze the results of Wang et al. (2009) to show that SN la with different 
Si II /i6355 absorption velocities at maximum light have different intrinsic colors. 
Accounting for these intrinsic color differences reduces the scatter of Hubble 
residuals by ~ 30%, while using SN la from a "normal" subsample reduces the 
scatter by ~ 40%. Although Foley & Kasen (2010) do not cross-validate their 
results or comment on their statistical significance, their analysis suggests that 
spectroscopy could be used to select a subsample of "well-behaved" SN la for 
more precise distance measurements. 
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Appendix A: Sampling variance of weighted mean 
square error 

If we assume the prediction errors are distributed normally with 
total variance, Eq. 5: A/i , ~ A^(0, cr^), then the sampling variance 
of the weighted mean square error is 

2 

Var[WRMS^] = j Yj ^^^^t ^^-^^ 



The standard error on WRMS is then VVarLWRMS2]/(2 x 
WRMS). 
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Appendix B: IVIaximum likeliliood estimators for 
intrinsic prediction error and covariance 

Let the intrinsic variance of predictions (cr^^^^) from method P, Q 
be cTp, cTg. Denote the intrinsic covariance between distance pre- 
diction errors from P and Q as cpq. The intrinsic correlation is 
ppQ = cpQ/io-pCTQ). Let = {cTp, ctq^ppq) be the vector of this 
triplet of quantities. These parameters can be arranged in an in- 
trinsic covariance matrix: 



I^intCf) = ( 

\PPQ0-P(^Q 



pPQcrpcrQ\ 
4 ) 



(B.l) 



We wish to estimate them from the cross-validated predictions. 
We derive estimators for this intrinsic covariance using maxi- 
mum UkeUhood, and also derive its uncertainty. 

We assume that the pair of prediction errors Aju^ = 
(A/if , A/<^) are jointly distributed normally around zero with to- 
tal covariance 



i:,(0) = i:; + i:;ec + ^int(0). 



(B.2) 
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The measurement error covariance matrix, S^, contains the mea- 
surement variances, cr^^, for model P and model Q, on the di- 
agonal, and any covariance due to observational error in the ofF- 
diagonal. Since supernova s is subject to the same random pe- 
culiar velocity under both models P and Q, the covariance from 
pecuUar velocity dispersion is 




(B.3) 



Appendix C: Results for other spectroscopic 
indicators at maximum liglit 

We present our results using the absorption velocity (Vabs). the 
full-width at half-maximum (FWHM), the relative absorption 

depth (dabs), the pseudo-equivalent width (pEW), and the vari- 
ous spectroscopic ratios [Eqs. 16-22] in Tables C.1-C.5. 



The negative log likelihood for the unknown (crp, o-q,ppq), 
given the set of distance predictions is 



1 

-m = 2 2 log det \_2ni:Aff)\ + A//[5:,(0)A//, (B.4) 

We numerically maximize the UkeUhood with the constraints 
(Tp, ctq > and \ppq\ < 1. Once we have found the maximum 
likelihood estimate (MLE) 0, we can compute its error by nu- 
merically evaluating the Hessian of the negative log hkehhood, 
H{ff). The sampling covariance (error) of the MLE is estimated 
from the Fisher information: = H^^(ff). The standard errors 
in each of (crp, o-q,ppq) are the square roots of the diagonal el- 
ements of Fg. The off-diagonal elements contain the estimation 
covariance between the three parameters. If the difference in in- 
trinsic prediction error between the two models is A = crp - erg, 
the sampUng variance of A is 

Var[A] = Var[o-p] - 2Cov[crp, ctq] + Var[crg] 

e e e 

where F^^' is the (;, /') element of the error covariance matrix 
of the MLE. This error estimate accounts for covariance from 
random peculiar velocities and the intrinsic correlation between 
two models. Notably, a large \ppq\ will affect the significance of 
the difl'erence, A. 

From the prediction errors of a single method, {A// J, we can 
estimate the rms intrinsic prediction error o-pred- The negative log 
UkeUhood simplifies to 



1 ^ Aju^ 

(B.6) 

The maximum likelihood estimate cr^^^^ is found by minimizing 
this or finding the zero of the score function L'icr^^^) = 0. If 
N is large enough, the standard error on d^^^ can be estimated 
using the Fisher information at the MLE: 



-^"(O^red) - Z .^2 + ^2 ^ 



^ 1 



(B.7) 

An estimate of the sampling variance of the maximum UkeUhood 
estimate of the intrinsic variance is the inverse of the Fisher in- 
formation Var((T2 ^j) - [-L"((Tpj.gj)]"'. The standard error of 

o'pred itself is the square root of Var[a-p^gj]/(4a-pj.g^). This esti- 
mate of the intrinsic dispersion "subtracts" out the contribution 
of random peculiar velocities and measurement error to the total 
dispersion. 



21 



S. Blondin et al.: Do spectra improve distance measurements of Type la supemovae? 



Table C.l. Vabs (units of l(f km s at maximum light from 10-fold CV 



Line 


r 


WRMS 


C^pred 








Vabs 


SilU4130 


-0.19 + 0.20 


0.289 ± 0.035 


0.271 ±0.038 


0.59 ±0.12 


0.073 ± 0.037 (2.0cr) 


33 


S1U5454 


0.33+0.16 


0.261 ± 0.037 


0.241 ± 0.040 


0.65 ±0.13 


0.061 ±0.038 (1.6cr) 


25 


Si 11/15972 


-0.02 ±0.19 


0.301 ± 0.040 


0.282 ± 0.043 


0.47 ±0.17 


0.106 ±0.048 (2.20-) 


28 


SilU6355 


-0.60 ±0.13 


0.275 ± 0.033 


0.256 ± 0.035 


0.73 ± 0.09 


0.090 ±0.031 (2.90-) 


35 


ta.Vabs) 


SiIU4130 


-0.11 ±0.20 


0.297 ± 0.037 


0.281 ±0.039 


0.68 ±0.10 


0.081 ±0.035 (2.3a-) 


33 


SIU5454 


0.20 ± 0.21 


0.257 ± 0.036 


0.234 ± 0.040 


0.71 ±0.12 


0.054 ±0.034(1.60-) 


25 


SilU5972 


0.44 ± 0.25 


0.305 ± 0.041 


0.288 ± 0.043 


0.50 ±0.16 


0.107 ±0.046 (2.30-) 


28 


SilU6355 


-0.41 ± 0.08 


0.261 ±0.031 


0.242 ± 0.034 


0.75 ± 0.09 


0.073 ± 0.028 (2.60-) 


35 


(C, Vabs) 


Si 1U4130 


0.78 ±0.14 


0.229 ± 0.029 


0.210 ±0.031 


0.71 ±0.10 


0.008 ± 0.029 (0.3cr) 


33 


SIU5454 


0.75 ± 0.22 


0.240 ± 0.035 


0.222 ± 0.037 


0.82 ± 0.08 


0.038 ± 0.025 (1.50-) 


25 


Sini5972 


0.20 ±0.10 


0.237 ± 0.032 


0.217 ± 0.035 


0.69 ±0.11 


0.035 ±0.030(1.20-) 


28 


SilU6355 


0.56 ±0.13 


0.227 ± 0.028 


0.207 ± 0.030 


0.70 ±0.10 


0.035 ±0.026(1.30-) 


35 


(X1,C, Vabs) 


SilU4130 


0.90 ±0.11 


0.209 ± 0.026 


0.188 ±0.029 


0.90 ± 0.03 


-0.013 ±0.017 (0.8tT) 


33 


S1U5454 


0.62 ± 0.24 


0.217 ±0.031 


0.193 ±0.034 


0.96 ± 0.02 


0.012 ± 0.012 (l.Otr) 


25 


Si 1U5972 


0.72 ±0.15 


0.204 ± 0.027 


0.179 ±0.031 


0.87 ± 0.05 


-0.003 ±0.021 (0.1 cr) 


28 


SilU6355 


1.07 ±0.12 


0.179 ±0.021 


0.151 ±0.025 


0.83 ± 0.06 


-0.020 ±0.019(1.10-) 


35 



Table C.2. FWHM (units of 10^ A) at maximum light from 10-fold CV 



Line 


r 




WRMS 




Pxi,c 










FWHM 


SiIU4130 - 


-1.28 ± 


0.09 


0.297 


± 


0.036 


0.277 


± 


0.040 


0.57 


± 


0.14 


0.093 


± 


0.039 (2.40-) 


32 


S1U5454 


0.15 ± 


0.57 


0.957 


± 


0.126 


0.952 


± 


0.126 


-0.01 


± 


0.22 


0.769 


± 


0.129 (6.O0-) 


29 


SilU5972 


0.14 ± 


0.20 


0.434 


± 


0.059 


0.424 


± 


0.059 


0.47 


± 


0.15 


0.246 


± 


0.063 (3.90-) 


28 


SilU6355 - 


-0.32 ± 


0.07 


0.279 


± 


0.033 


0.259 


± 


0.036 


0.72 


± 


0.09 


0.089 


± 


0.031 (2.90-) 


35 


(xi,FWHM) 


Sini4130 - 


-1.34 ± 


0.09 


0.293 


± 


0.035 


0.273 


± 


0.039 


0.49 


± 


0.15 


0.094 


± 


0.037 (2.50-) 


32 


Sn^5454 


0.44 ± 


1.17 


0.715 


± 


0.092 


0.708 


± 


0.093 


0.13 


± 


0.21 


0.524 


± 


0.084 (6.20-) 


29 


SiIU5972 


0.33 ± 


0.24 


0.507 


± 


0.068 


0.499 


± 


0.069 


0.40 


± 


0.15 


0.318 


± 


0.073 (4.40-) 


28 


SiIU6355 - 


-0.24 ± 


0.08 


0.276 


± 


0.033 


0.256 


± 


0.035 


0.75 


± 


0.09 


0.085 


± 


0.029 (2.90-) 


35 


(cFWHM) 


SiIU4130 - 


-0.89 ± 


0.14 


0.258 


± 


0.032 


0.236 


± 


0.035 


0.73 


± 


0.10 


0.052 


± 


0.031 (1.70-) 


32 


S1U5454 


0.50 ± 


0.44 


0.495 


± 


0.065 


0.485 


± 


0.066 


0.22 


± 


0.21 


0.304 


± 


0.073 (4.2cr) 


29 


SiIU5972 


0.28 ± 


0.11 


0.359 


± 


0.049 


0.348 


± 


0.050 


0.60 


± 


0.14 


0.167 


± 


0.050 (3.3cr) 


28 


SiIU6355 - 


-0.04 ± 


0.09 


0.247 


± 


0.030 


0.227 


± 


0.032 


0.79 


± 


0.07 


0.054 


± 


0.025 (2.20-) 


35 


(xi,c,FWHM) 


SiIU4130 


0.24 ± 


0.51 


0.233 


± 


0.030 


0.214 


± 


0.031 


0.93 


± 


0.03 


0.031 


± 


0.017 (1.8a-) 


32 


SIU5454 


0.01 ± 


0.41 


0.331 


± 


0.043 


0.315 


± 


0.046 


0.43 


± 


0.17 


0.135 


± 


0.047 (2.9cr) 


29 


SiIU5972 


0.26 ± 


0.08 


0.362 


± 


0.050 


0.352 


± 


0.050 


0.66 


± 


0.10 


0.170 


± 


0.045 (3.80-) 


28 


SiIU6355 


0.28 ± 


0.03 


0.193 


± 


0.023 


0.168 


± 


0.026 


0.90 


± 


0.04 


-0.006 


± 


0.015 (0.4a-) 


35 
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Table C.3. rfabs at maximum light from 10-fold CV 



Line 


r 


WRMS 


C^pred 










SilU4130 


-2.51+0.29 


0.283 + 0.035 


0.267 + 0.037 


0.54 ±0.11 


0.076 ±0.036 (2.1cr) 


33 


S1U5454 


-3.07 + 0.74 


0.280 + 0.039 


0.259 + 0.043 


0.37 ± 0.20 


0.063 ±0.050 (I.3cr) 


25 


SilU5972 


-3.54 ± 1.31 


0.357 ± 0.046 


0.336 ±0.051 


0.34 ± 0.20 


0.160 ±0.056 (2.90-) 


28 


SilU6355 


-1.66 + 0.18 


0.277 + 0.033 


0.257 + 0.036 


0.57 + 0.13 


0.088 + 0.032 (2.70-) 


35 


(Xi.dabs) 


SiIU4130 


-3.84 ± 1.60 


0.384 + 0.048 


0.375 + 0.048 


0.24 + 0.18 


0.179 + 0.052 (3.40-) 


33 


SIU5454 


-2.53 + 0.65 


0.270 ± 0.037 


0.248 ± 0.042 


0.52 ±0.17 


0.052 ±0.043(1.20-) 


25 


SiIU5972 


-3.13 ± 1.12 


0.260 ± 0.038 


0.239 ± 0.042 


0.34 ± 0.22 


0.023 ±0.022(1.00") 


28 


SilU6355 


-0.77 ± 1.14 


0.306 ± 0.037 


0.290 + 0.039 


0.62 + 0.12 


0.119 + 0.037 (3.20-) 


35 


(C,4bs) 


Si 1U4130 


-2.06 + 0.27 


0.230 + 0.029 


0.210 ±0.031 


0.80 ± 0.08 


0.019 ± 0.026 (0.7cr) 


33 


Sn/15454 


-3.90 ± 0.33 


0.228 ± 0.032 


0.204 ± 0.036 


0.49 ±0.17 


0.008 ± 0.038 (0.2cr) 


25 


SiIU5972 


-2.65 + 0.70 


0.243 + 0.031 


0.214 + 0.037 


0.60 + 0.13 


0.044 + 0.037(1.20-) 


28 


SiIU6355 


-1.43 + 0.14 


0.226 + 0.027 


0.203 + 0.030 


0.78 + 0.08 


0.031 +0.023(1.30-) 


35 


(Xl,C,rfabs) 


SilU4130 


-1.87 + 1.04 


0.248 +0.031 


0.231 ±0.033 


0.87 ± 0.05 


0.039 ± 0.020 (2.0cr) 


33 


S1U5454 


-3.26 + 0.35 


0.203 + 0.028 


0.175 ±0.033 


0.75 ±0.10 


-0.022 ± 0.030 (0.7tr) 


25 


Si 11/15972 


-2.24 + 0.55 


0.217 + 0.028 


0.188 ±0.033 


0.69 ±0.12 


0.014 ±0.030 (0.5cr) 


28 


SilU6355 


0.18 ±0.67 


0.219 ± 0.026 


0.198 ±0.029 


0.97 ± 0.01 


0.028 ±0.011 (2.50-) 


35 



Table C.4. pEW (units of 10^ A) at maximum light from 10-fold CV 



Line 


r 




WRMS 
















pEW 


SiIU4130 


-2.09 ±0.18 


0.257 


± 


0.032 


0.238 


± 0.034 


0.65 


± 


0.12 


0.041 


± 


0.033 (1.20-) 


33 


FelU4300 


-0.20 H 


tO.23 


0.298 


± 


0.035 


0.282 


± 0.037 


0.64 


± 


0.09 


0.107 


± 


0.035 (3.10-) 


36 


FelU4800 


-0.32 H 


t0.05 


0.269 


± 


0.032 


0.250 


± 0.034 


0.67 


± 


0.11 


0.076 


± 


0.031 (2.5cr) 


36 


S lU.i5454,5640 


-0.65 H 


tO.26 


0.269 


± 


0.035 


0.249 


± 0.038 


0.63 


± 


0.13 


0.076 


± 


0.037 (2.1fr) 


29 


SiIU6355 


-0.55 H 


t0.05 


0.256 


± 


0.030 


0.235 


± 0.033 


0.73 


± 


0.09 


0.058 


± 


0.029 (2.0O-) 


35 


(xi,pEW) 


Sini4130 


-2.53 ± 0.28 


0.262 


+ 


0.032 


0.244 


+ 0.035 


0.63 


+ 


0.12 


0.046 


+ 


0.033 (1.40-) 


33 


Feni4300 


-0.13 ±0.26 


0.297 


± 


0.035 


0.281 


± 0.037 


0.69 


± 


0.10 


0.104 


± 


0.034 (3.10-) 


36 


FeIU4800 


-0.37 H 


t0.05 


0.258 


± 


0.030 


0.239 


± 0.033 


0.78 


± 


0.08 


0.064 


± 


0.026 (2.50-) 


36 


S lU.i5454,5640 


-0.16 H 


tO.62 


0.288 


± 


0.038 


0.272 


± 0.040 


0.68 


± 


0.12 


0.096 


± 


0.040 (2.40-) 


29 


SilU6355 


-0.43 ± 0.06 


0.265 


+ 


0.031 


0.245 


+ 0.034 


0.74 


+ 


0.09 


0.067 


+ 


0.029 (2.30-) 


35 


(cpEW) 


SilU4130 


-1.57 H 


tO.24 


0.225 


± 


0.028 


0.205 


± 0.030 


0.94 


± 


0.03 


0.006 


± 


0.014 (0.4O-) 


33 


FeIU4300 


0.14 H 


t0.17 


0.240 


± 


0.029 


0.221 


±0.031 


0.84 


± 


0.05 


0.043 


± 


0.021 (2.0O-) 


36 


FeIU4800 


0.23 ± 0.08 


0.230 


+ 


0.028 


0.212 


+ 0.029 


0.81 


+ 


0.06 


0.035 


+ 


0.022 (I.60-) 


36 


SIUA5454,5640 


-0.90 ±0.14 


0.223 


+ 


0.030 


0.201 


+ 0.032 


0.81 


+ 


0.08 


0.026 


+ 


0.024(1.10-) 


29 


SiIU6355 


-0.31 ± 0.05 


0.231 


± 


0.028 


0.210 


± 0.030 


0.89 


± 


0.04 


0.029 


± 


0.018(1.60-) 


35 


(xi,c-,pEW) 


SilU4130 


0.01 H 


tO.78 


0.249 


± 


0.031 


0.232 


± 0.033 


0.99 


± 


0.00 


0.022 


± 


0.008 (2.70-) 


33 


FelU4300 


0.25 H 


t0.16 


0.216 


± 


0.026 


0.196 


± 0.028 


0.99 


± 


0.00 


0.013 


± 


0.007 (1.90-) 


36 


FelU4800 


0.17 H 


t0.09 


0.212 


± 


0.026 


0.193 


± 0.027 


0.97 


± 


0.01 


0.013 


± 


0.010(1.30-) 


36 


S IU.i5454,5640 


-0.32 H 


tO.22 


0.210 


± 


0.028 


0.186 


±0.031 


0.99 


± 


0.01 


0.013 


± 


0.009(1.40-) 


29 


SiIU6355 


0.31 ±0.13 


0.215 


± 


0.026 


0.193 


± 0.028 


0.95 


± 


0.02 


0.013 


± 


0.011 (1.20-) 


35 



S. Blondin et al.: Do spectra improve distance measurements of Type la supemovae? 
Table C.5. !??(Ca) and KiSi) at maximum light from 10-foId CV 



Line y WRMS o>ed Pxi^ A^i£ Nsm^ 

;Rx 

«(CaS) -0.10 + 0.03 0.288 ±0.043 0.273 +0.046 0.52 ±0.15 0.085 ± 0.046 (1.8cr) 22 

«(Si) -2.30 ±0.45 0.282 ±0.037 0.259 ±0.041 0.31 ±0.20 0.097 ± 0.048 (2.0cr) 28 

«(SiS) 1.52 ±0.13 0.286 ±0.038 0.270 ± 0.040 0.69 ±0.11 0.1 1 1 ± 0.040 (2.8cr) 29 

■K(SiSS) 1.24 ±0.18 0.260 ±0.031 0.242 ± 0.033 0.73 ± 0.09 0.053 ± 0.029 (I.80-) 36 

??(S,Sif 0.09 ±0.06 0.263 ± 0.036 0.242 ± 0.040 0.59 ±0.15 0.065 ± 0.039 (1.7cr) 26 

!/?(Si,Fe) -1.69 ±0.23 0.275 ± 0.036 0.257 ±0.038 0.53+0.15 0.079 ± 0.041 (1.9o-) 29 

{xuHx) 

•K(CaS) -0.15 ±0.03 0.273 ±0.041 0.259 ± 0.043 0.63 ± 0.14 0.071 ± 0.042 (1.7o-) 22 

«(Si) -2.39 ±0.20 0.196 ±0.033 0.175 ±0.035 0.27 ± 0.24 0.026 ± 0.045 (0.5cr) 28 

«(SiS) 1.56 ± 1.52 0.398 ±0.052 0.386 ± 0.053 0.35 ±0.17 0.214 ± 0.067 (3. 2fr) 29 

«(SiSS) 1.34 ±0.33 0.277 ± 0.033 0.260 ± 0.035 0.70 ± 0.09 0.071 ± 0.029 (2.4fr) 36 

«(S,Sif -0.28 ±0.13 0.278 ± 0.038 0.257 ± 0.041 0.64 ± 0.12 0.079 ± 0.036 (2.20-) 26 

!y?(Si,Fe) 0.09 + 0.85 0.285 + 0.038 0.269 + 0.040 0.64 + 0.12 0.091 + 0.038 (2.4o-) 29 

iCRx) 

•R(CaS) 0.08 ±0.03 0.248 ±0.037 0.234 ±0.040 0.73 ±0.11 0.043 ± 0.035 (1.2o-) 22 

■«(Si) -1.68 ±0.19 0.190 ±0.025 0.158 ±0.030 0.59 ±0.13 -0.007 ± 0.030 (0.2tr) 28 

7?(SiS) 0.92 ±0.21 0.240 ±0.032 0.220 ± 0.034 0.86 ± 0.06 0.062 ± 0.022 (2.8cr) 29 

■R(SiSS) 0.81 ±0.19 0.229 ±0.027 0.208 ± 0.030 0.95 ± 0.02 0.021 ± 0.01 1 (1.9cr) 36 

■^(8,81)" 0.25 ±0.04 0.220 ±0.032 0.200 ± 0.033 0.88 ± 0.05 0.017 ± 0.022 (0.80-) 26 

■^(Si,Fe) -1.99 + 0.35 0.216 + 0.029 0.196 + 0.031 0.88 + 0.05 0.018 + 0.021 (0.9o-) 29 

{xuc^Hx) 

!/?(CaS) 0.04 + 0.02 0.215+0.033 0.202 ±0.035 0.99 + 0.01 0.006 + 0.012 (0.5a-) 22 

7?(Si) -1.23 + 0.87 0.209 + 0.027 0.179 + 0.032 0.79 + 0.08 0.021 + 0.026 (O.80-) 28 

«(Si8) -0.98 ±0.20 0.205 ±0.028 0.184 ± 0.030 0.84 ± 0.06 0.022 ± 0.017 (1.3o-) 29 

■/?(Si88) -1.06 ±0.29 0.225 ± 0.027 0.205 ± 0.029 0.96 ± 0.02 0.013 ± 0.013 (l.Otr) 36 

'R(S,8i)'' -0.11 ±0.09 0.216 ±0.031 0.195 ±0.033 0.99 ± 0.01 0.013 ± 0.011 (1.2o-) 26 

!y;(Si,Fe) -0.16 + 0.90 0.233+0.032 0.216 + 0.033 1.00 + 0.00 0.035 + 0.009 (3.9o-) 29 



Notes. <°>Infact??(S,Si)/10. 
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